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1.0 SUMMARY 


Significant differences exist in the noise generated by engines in flight and 
engines operating on the test stand. It has been observed that these differ- 
ences are reduced by the use of an inflow control structure in the static test 
situation. It is a purpose of this contract (NASU15085) to produce a design 
systen for inflow control structures. This report first describes the results 
of a test program conducted to develop various theoretical models used in the 
design system. Following this the design system is developed. Finally an as- 
sessment of this procedure Is made using measured data. 

The primary result of this work is a step by step Inflow control structure de- 
sign procedure in terms, for the most part, of parameters quantifiable by the 
designer. The assessment of this procedure using blade mounted transducer data 
was inconclusive due to the transducer responding to stimuli other than the 
engine inflow field. However an assessment made with hot film data was very 
encouraging. 


2 0 INTRODUCTION 


Reduction of Jet engine noise levels Is a continuing evolutionary ^ocess and 
significant noise reduction has been achieved since the advent of the original 
cofiwerclal jet aircraft in the early. I950's. Engine noise reduction features 
are usually verified by means of static or flight testing. Since the cost of 
flight testing Is much greater than that of static testing, use of the latter 
technique results in significantly lower costs being accrued to the develop- 
ment of engine noise reduction features. As such, the use of static testing 
techniques should be fully exploited. The purpose of this contract Is to de- 
velop, for turbofan engines, static testing techniques that result in measure- 
ments of fan noise that are representative of Inflight levels. 

It has been noted by various observers that the noise produced by the fan of 
turbofan engines operating statically on the test stand is greater than that 
produced when the engine Is operating under flight conditions. This may be a 
result of both fan blade passing tone and broadband levels being contaminated 
by extraneous noise sources present during static testing. As a consequence, 
predictions of flight noise levels using static data are consistently high. A 
no less Important ramification is that flight noise sources cannot be identi- 
fied from noise measurements from statically operating engines. Correspond- 
ingly, noise reduction techniques cannot be evaluated on the test stand tn the 
presence of a contaminating noise floor. To identify the source of extraneous 
noise it is necessary to note that, stvticatly, the fan inteiacts with a more 
distorted Inflow field than it does when In flight. There are several distinct 
characteristics of the Inflow field that are quite different when the engine 
is operating statically, each of which could produce extraneous noise. 

Firstly, the ambient atmospheric turbulence field in the vicinity of the 
ground is more turbulent than at higher altitudes. In addition, this turbul- 
ence field Is convected through a very high flow contraction whefi the engine 
Is operating statically, whereas. In flight, the turbulence field convects 
through a very small contraction on its way to the fan. This high flow 
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contraction In tho static case resolts In a distortion of the field In which 
the "eddies" are "stretched", then "chopped" by successive fan blades, produc- 
Z “buTsts" of discrete noise that are virtually absent In "he flight opera- 

tion of the engine, 

secondly, 'n the static case, there esists stand structure and a ground plane 
that are sources of flov disturbances. Engine Ingested air passes ov r t 
stand structure and the ground resulting In vortices and " 

field. Engine case protuberances can also generate Ingestib e 
Usually these sources of distortion do not e«1st Inflight since the air 
gested into an engine has not passed over any boundaries. 

Thirdly, the nacelle boundary layer Is different during static ^ 

the engine. In flight, at takeoff and approach 

r^girlhan re :an. the locus of the stagnation to e -r the 

leaLg edge of the nacelle. In static operation, however, air Is Ingeste 
ir a wuerwider strea^tube causing the stagnation locus to move rear*, d on 
the outside of the flight nacelle. The acceleration of the flow 
stagnation locus around the lip region produces higher local ‘ " 

those occuring In flight and consequently, a larger adverse 
along the nacelle Inner wall than that which Is encountered In 
boundary layer then becomes thicker In this region and may even separate (tbe 
use oJ belTmouth Inlets during static testing Is an attempt to ai.evi h v 
situation). Due to these differences In boundary layers, the area of J 
t o that affected by the boundary layer during static operation Is probab y 
t 1: :.cene boLary layer In static operation Is affected y the 
presence of free stream flew distortions which can locally PJturb t. 

Lans probable, then, that the nacelle boundary layer generated by 
operation of the engine constitutes a larger flow disturbance than Its flight 
counterpart and thus could be a cause of higher noise levels. 
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Finally, since it is possible for the flow to be drawn from all angles by an 
engine on the test stand, turbulent flows from the jet plume may be reingest- 
ed. This distortion source is not present in flight. These disturbancos~(i.e., 
atmospheric turbulence, ground plane and stand Induced distortions, dissimUar 
nacelle boundary layer and jet plume reingestion) are considered to be the 
most important extraneous noise sources in static engine operation. In order 
to Obtain useful static acoustic data, it is- therefore necessary to develop 
techniques which modify the inflow field so that, for all. intents and pur- 
poses, the fan is operating as it would in flight. In the past, various tech- 
niques for accomplishing this simulation have been used. Including mounting 
engines in wind tunneU and using devices upstream of the fan to condition the 
inflow. Inflow Control Structures (ICS) for conditioning the flow have been 
mounted upstreaiii of the engine by several investigators. A boundary layer suc- 
tion system has also been used In conjunction with a screen. Both of these 
techniques have resulted in reduced radiated noise levels. Indicating the re- 
duction of inflow distortion. 

In view of the encouraging results achieved by the use of ICS's, the present 
contract was awarded for the purpose of developing an inflow control screen 
design procedure and a flight noise prediction procedure using data gathered 
from the static testing of engines equipped with such a structure. The use of 
wind tunnels, or other techniques for simulating the flight environment are 
not considered under this contract. In addition, the problem of nacelle bound- 
ary layer modification is not addressed under this contract, but is certainly 
, . an important area for future work. The major components of the contract are: 

Phase I Definition of Atmospheric Turbulence Characteristics and Engine 
Sensitivity Study. 


Phase II Development of Inflow Con<^*'oI Structure (ICS) Preliminary Design 
Systan. 


Phase III 


Interim Procedures Report Development and Coordination. 
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This report describes e prelhninery design system for the Inflw control- 
structure (ICS). Specifically Tasks F, G and H are addressed here. 

An interim Phase II report, previously published, contains details of the ana- 
lytical models used In the design system. The results of a test pror«" de- 
sired to assess and/or modify these analytical models are reported here. Sub- 
seguently, the models arc modified and used to develop an ICS design system. 
This desi systen Is exercised for the JTgO and JT150. Finally, the ICS de- 
sign system was assessed using blade mounted pressure transducer-(BMT) data 
gathered under the P&WA/Boetng Joint Noise Reduction Program. 

This report, in conjunction with the Interim Phase II report conprises the 
final Phase II report. It should be noted, that in this report and the Interim 
Phase II report, the term "screen" is used in ifs general sense, thus, 
perforated plate, honeycomb and gauze are all types of screen. 


3.0 SYMBOLS 


a 

Sq 

B 

b 

D 

Do- 

d 

F 

F 

f 

AH 

K 

KE 

k 

L 

C 

M 

N 

NR 

P 

P 


Characteristic Dimension 
Speed of Soun^ 

Fan Blade Nunber 
Wake Width 
Contraction Length 
Inlet Ditineter 
Detail Dimension 
Velocity PSD 
Normalized Velocity PSD 
Frequency 

Total Pressure Loss 
Pressure Drop Coefficient 
Kinetic Energy 
Wavenumber Vector 
Integral Length Scale 
Thickness 
Contraction Ratio 
Mach Number 

Isotropic Turbulence Characteristic 
Noise Reduction 
Pressure PSD 
Pressure 

Function Defined in 6.4 


q Dynamic Head 


R 

Re 

r 

SPL 

s 

U 

AUi 

u 

X 

Xo 

a 

2 


zq 

2REF 

a 

P 

r 

r: 

y 

y* 

e 


M 


V 

p 

a 

0 

w 


Radius 

Reynolds Number 
Vortex Radius 
Sound Pressure Level 
Defined In 6.4 
Mean Velocity Vector 
Velocity Deficit 
Turbulent Velocity Vector 
Streamwise Distance 
Virtual Source Distance 
Position Vector 
Engine Height 

Terral n Xbaracteri sti c Hel ght 
Reference Height 
Flow Angle Ratio 
Bias 

Circulation 

Velocity Spectrum Tensor 
Isotropic Turbulence Characteristic 
Blade Relative Mean Flow Angle 
Angie 

Minimum Sensitive Scale 

Mean Square Velocity Component Ratio 

Kinematic Viscosity 

Air Density 

Open Area Ratio 

BMT Path Angle 

Circular Frequency 
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Subscripts 
A 

AMB 
B 
c 

EFF 
F 

H 
I 

ICS 

i, j 

L 
m 

max 
NS 
n 
0 
P 
R 

REF 
S 

SC 
TU 
TD 

i 

1 


Pre, Contraction Or Screen 
Atmosphere 

Ambient 

Post, Contraction or Screen 
Core 

Effecti ve 
Flight 

Final Contraction 
Honeycomb 

Initial Contraction 

At ICS 

Indices 

Lower 

Microphone 

Maxlmim 

No Screen 

Upwash Component 

Fan Face Conditions 

Perforated Plate (Gauze) 

Recel ver 

Reference 

Source, Static 

Screen 

Total Upstream 
Total Downstream 
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4.0. .BACKSROIMD 


It has long been known that the flight and static test envlroments of an en- 
gine are different. For Instance. In 1966, Sofrin and McCann (I) recognized 
the existence of noise generation mechanisms peculiar to ground operation of 
an engine. They Identified "natural atmospheric air currents" and ground vor- 
tices as causes of Inflow distortion that, upon Interaction with a rotor, pro- 
duce noise. They also surmised that the Inflow Into a fan In flight was not 
distortion free, but probably had a different distortion flow field than which 
existed during static testing. 

Subsequently, accumulated evidence {2. 3. 4. 5. 6) showed that not only were 
static and.f light inflow distortions different, but that the static distort on 
field generated noise levels, especially blade passing tone levels (and possi- 
bly Its harmonics) that were significantly greater than those generated in 
flight. However, while the source of the extraneous static noise Is generally 
ascribed to the more disturbed Inflow field, the particular distortion type 
responsible Is not clearly defined. For example. Cunpsty and Lowrie (2) pro- 
duced results that indl-cated that boundary layer changes were Important at re- 
lative tip Mach numbers less than .85. In the flow outside the boundary l^yer. 
they suggest that distortions produced by the contraction of the ambient at- _ 
mospheric turbulence field are Important. Hanson’s (7. 8. 9) measurements of 
Inlet flow fields are explained by that author also as the result of contrac- 
tion of the atmospheric turbulence, although, in some situations, he also 
Identifies distortions due to flow over the stand structure. Distortions have 
also been detected that were generated by protuberances on the exterior 
nacelle casing (10). Earlier work by Povinelll et al (II) Identified distor- 
tions of the Inflow by the installation structure. Hodder (12) examined among 
other things, the noise generated by a JT15D at two heights from the ground, 
and found the ground vortex Interacting with the rotor to be a significant 
source. Recent P&WA/Boeing experience under the Joint Noise Reduction Program 
(JNRP) indicates that on a full size test stand, ground effects, stand struc- 
ture and possibly abnc-pherlc turbulence contribute to inflow distortions 


(13). The above review indicates that any or all of the distortion sources can 
exist but to a varying degree depending on the details of the test stand con- 
figuration. 

Significant suppression of these flow distributions has been achieved by the 
mounting of Inflow Control Structures (ICS) upstream of the inlet. These 
structures consist of screening materiaJs either supported by a frame or sub- 
stantially self supporting with the minimum^of unnecessary structure. 

Published- experience with ICSs ivas been quite varied. The scales of the fans 
has ranged from model to the large high bypass ratio fan. The ICS has been 
placed in the inlet duct and externally. 

The in-duct ICS has consisted of either a honeycomb panel mounted upstream of 
the rotor (12* 14, 15) or multiple gauze screens at the inlet entrance (4). 

All except Leggat and Siddon observed a reduction in tone level, however, 
these investigators did note a change in blade loading caused by the suppres- 
sioa of transverse velocities in a ground vortex. Suppression of this vortex 
nearer its stagnation point resulted in noise reductions. In conjunction with 
the noise reduction, Hodder (12) observed a reduction in axial length scale, 
as determined by making auto correlations of the hot wire signal obtained at 
the fan face. Hodder ascribed the noise reduction to a change in length scale 
caused by the screen. However, Jones, et al, (15) demonstrated that honeycomb 
screens mounted In the inlet duct have an acoustic transmtsssion loss associ- 
ated with then, so there is some question as to the reason for the noise re- 
ductions when the screens are mounted in the inlet. 

The externally mounted ICS's do not suffer from this problem of acoustic 
transmission loss to any great degree (15, 16). In general, all screens of 
this type (17. 18. 15, 19, 20, 21. 13, 22, 23, 24) produced a significant re- 
duction in tone noise at subsonic relative tip Mach numbers, whether in the 
anechoic chanber. wind tunnel or on the test stand. The work of Rogers ot a I 
(13), Atvars et al (22) and McArdle et al (23) demonstrates conclusively, with 
the use of blade mounted pressure transducers, the substantial clean up of the 
inflow field that can be achieved with the ICS. In addition, Rogers and Ganz 


in 






(13). by compiirliig st.tic and fll^t BUT data show that the ICS goes a long 
way to reproducing the flight fan face environment. Woodward et al (24), and 
Kantola et al (21). have also used inflow boundary layer modification in con- 
junction with an Inflow Control Structure to produce encouraging results. 

While universal screen design criteria have not yet been determined, certain 
desirable characteristics have emerged based on the literature reviewed. These 
Characteristics are sumnarized below; 

0 The minimization of acoustic transmission loss is enhanced by 
use of. external ICSs as opposed to use of ICSs mounted in the 

engine inlet (15, 16). 

0 Distortion suppression, primarily of transverse velocities, is 
effected by use of honeycomb. Honeycomb characteristics sug- 
gested by Ho, et al, (25) areJthat a length to diameter ratio 
of 5 to 10 is necessary. However, Ginder (26) shows that in- 
creasing this ratio beyond 2 does not improve disturbance sup- 
pression. This latter view is born out, in seme degree, by 
Prandtl (27). 

0 Axial velocity distortion suppression is effected by use of gauzes or 
perforated plate because of their relatively high pressure loss coef- 
ficients. 

0 Distortions generated by honeycomb structure can be suppressed by use 
of gauzes mounted downstream of the honeycomb, as suggested by the 
work in references (17), (25), (26) and (28). 

0 Gauzes and perforated plates can be mounted upstream of the honeycomb 
to protect it from damage by incoming debris. 



0 Inf low control devices should be designed to avoid large structural 
members and discontinuities (as noted by- Lowrie-and Newby (16) and 
PWA/Boeing experience) since they can give rise to self generated 
distortions entering the fan. As can be seen from the above review of 
the literature, the feasibility has been shown of using an Inflow 
control structure during static testing to effect a better simulation 
of the fan Inflight flow field. Certain desirable characteristics of 
ICS designs have been Identified, but as yet, a design procedure does 
not exist. An ICS design procedure that Is generally appplicable Is 
developed In the present work. 
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5.0 TEST PROGRAM (PHASE II TASK F) 


The purpose of this program was to provide data to empirically verify or cor- 
rect the analytical models for the ICS design procedure Identified In the In- 
terim Phase II report (29). The tests conducted were basically of two types, 
aerodynamic and acoustic. 


The aerodynamic tests consisted of determining the effect of flow contraction 
and screening on Inflow distortion elements and turbulence. In addition screen 
pressure drop coefficients were measured. 

The acoustic portion of the program consisted of measuring the effect of 
structural members and panel corners on the acoustic radiation field. (This 
was described In part In (29). Acoustic transmission loss measurements across 
various screens were also made. 

In the following subsections the aerodynamic measurements are discussed In 
Section 5.1 and the acoustic data In Section 5.2. 

5.1 The Effect of a Flow Contraction and Screening on Turbulence and Steady 
Flow Distortion 


The tests were carried out at Pratt & Whitney Aircraft's 25.4 cm suction rig 
facility. Flow through this rig is provided by three 450 hp vacuum pumps. Duct 
velocities ranging from 18.29 m/s to 121.95 m/s were generated to encompass 
the range of contraction ratios likely to be experienced by the engine on the 
test stand. A bellmouth was mounted on the suction duct to ensure a smooth 
transition of flow from the test room Into the 25.4 cm duct. A schanatlc dia- 
gram of the basic configuration Is shown In Figure 1. 

A turbulence generator was mounted upstream of the test room as shown In 
Figure 2, This consisted of a plenum which exhausted Into the test room 
through a series of opposing Jets (porcupine) (Figure 3), This arrang^ent 
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Figure J Turbulence Generating Device 


Screen 
Spec , _ 

Wake 

Gcn__ 

Turb. 

Gen 

Axial 

(OB).! 

Radial Positions 

-Cm)-, 

Hone 


Yes 

-20.32 

26.4 

55.88 

(CL to +10.16) 
(CL to *412.7) 


Duct Velocities Neas. 

jRanqe in m/s) L(Be_ 


(?4.39 to 121.95 
30.49 to 121.95 
30.49 to 121.95 


#2 

44 

#4 


IB Mono 


1C .953 cm H/C 
,953 cm H/C 
1.27 cm H/C 
1.27 cm H/C 
2.54 cm H/C 
2.54 cm H/C 
41X O.A. P.P. 
41X O.A. P.P. 
66» O.A. P.P. 
66X O.A. P.P. 

10 .953 cm 
Honeycomb 


1,27 cm 
Honeycomb 


2.54 cm 
Honeycomb 


4U O.A. 
Perf. Plate 


66t O.A. 
Perf. Plate 


3.81 Sq. Rod No 
5.08 X 7,62 Rod 
SYM. A/F (NACA 
Series 6) 

3.81 Sq, Hod 
3.81 Sq. Rod 
5.08 X 7.62 Rod 
5.08 X 7,62 Rod 
SYH. A/F (NACA 
Series 6) 

5.08 X 7.62 Rod 

- . Yes 


-20.32 

(-10.16 to +10.16) 

(30.49 to 121,96) 

#1 

55.88 

(-3.8 to +3.8) 

(30.49 to 121.95) 

f3 

25.4 

55,88 




25.4 

55.88 




25.4 



48.26 

(CL to +12.7) 

(30.49 to 121.96) 

#4 


30.48 

48.26 

30.48 

48.26 

30.48 

48.26 

30.48 

48.26 

30.48 


3.8 Sq. Rod No 

3.8 Sq. Rod 

5.08 X 7.62 Rod 
5.08 X 7.62 Rod 
Tip Vortex A/F 
(NACA 0012) 

Tip Vortex A/F 
(NACA 0012) 

3.81 Sq. Rod 
3.81 Sq. Rod 
5.08 X 7.62 Rod 
5.08 X 7.62 Rod 
Tip Vortex A/F 
(NACA 0012) 

Tip Vortex A/F 
(NACA 0012) 

3.01 Sq. Rod No 
3.81 Sq. Rod 
5.08 X 7.62 Rod 
5.08 X 7.62 Rod 
Tip Vortex A/F 
(NACA 0012) 

Tip Vortex A/F 
(NACA 0012) 

3.81 Sq. Rod 
3.81 Sq. Rod 
5.08 X 7.62 
Sq. Rod 

5.08 X 7.62 
Sq. Rod 

3.31 Sq. Rod No 
3.81 Sq. Rod 

5.08 X 7.62 Rod 

5.08 X 7.62 Rod 
Tip Vortex A/F 
(NACA 001?) 

Tip Vortex A/F 
(NACA 001?) 


48.26 (-12.7 to +12.7) 

30.48 

48.26 

30.48 

48.26 

30.48 

48.26 

30.48 

48.26 

30.48 

48.26 

30.48 

48.26 (-12.7 to +12.7) 

30.48 

48.26 

3C.48 

48.26 

30.48 

48.26 

30.48 

48.26 

30.48 

48.26 (-12.7 to +12.7) 

30.48 

48.26 

30.48 

48.26 

30.48 


(53.35 to 121.95) #3 


(53.35 to 121.95) #3 


(53.35 to 121.95) #3 


Figure 4 Configtirdtlons and Measurement 
Screening Tests 


Locations For Contraction and 


is based on a similar design of R. Betchov (30) and produces, to some degree, 
isotropic turbulence of reasonably high intensity. For the steady distortion 
part of this program the porcupine was removed and the outlet snout of the 
generator covered with honeycomb to provide a uniform flow over the distortion 
generating bodies. Measurements of the velocity fields werejnade 1n the duct 
itself and upstream of the suction duct. 

Outside of the duct, measurements were made with traversing hot wire probes. 
The allgiment of the probe traverse mechanism was accomplished by mounting the 
probes on support rails attached to 2.B4 cm thick steel bed plates. The probe 
positions were set using a laser beam directed along the rig axis. For the 
turbulence and vortex measurements, 'X' wires were used enabling two compon> 
ents of velocity to be determined. Specifically, these were the streamwise 
velxity component and the vertical transverse component. A single hot wire 
was used to measure the steady streamwise velocity distortions. 

Inside of the duct, in the high velocity region, the turbulence and mean velo- 
city measurements were made with 'X' hot wires and a wedge probe, respective- 
ly. The wedge probe allowed the magnitude and vertical angle of the component 
of the mean velocity vector in the vertical plane to be determined. Figure 4 
contains a table of the measurement types and probe locations for the contrac- 
tion and screening experiments. AXi turbulence fields were investigated with 
two probes to provide information about their transverse nature. A description 
of the measurement types and their accuracy is presented in Appendix A. 

5.1.1 The Contraction of Turbulence and Steady Distortions 

The contraction of turbulence was measured with the turbulence generating por- 
cupine in place, in the plenum. The measurements were made at three axial lo- 
cations through the contraction (55.8 an, 25.4 cm and -20.32 cm referenced to 
the beilmouth face with negative distance downstream) as shown in Figure 5. 

The two probes were separated by incremental amounts at each location in the 
horizontal transverse plane. The sum and difference of the signals from 
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Figure 5 Axial Measurement Locations for Turbulence Contraction Experiment 


each hotwire probe were recorded on magnetic tape and subsequently analyzed. 
The duct velocities examined In the contraction experiments nominally ranged 

frcm 30.48 m/s to 121.95 m/s. This allowed contraction ratios from about 1 to 
21.3 to be examined. 


Auto correlations, cross correlations, PSDs and mean square values were ob- 
tained directly fron the analysis. Figures 6 and 7 show how the auto correla- 
tions of the streainwlse and transverse velocity components change as the tur- 
bulence is convected through various contractions. It seems that the turbul- 
ence deviates somewhat from the Isotropic at the initial point since the auto 
correlation function of the streamwise component becmes negative; however, 
the turbulent kinetic energy is reasonably evenly distributed at this point. 
The streamwise component auto correlations at the two highest contraction 
ratios display sonewhat anomalous characteristics as the tendency to narrowing 
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Figure 





of the auto correlation with Increasing contraction is reversed in these two 
cases. This is explained later as due to transfer of turbulent kinetic energy 
fron the transverse velocity component by turbulent inertial effects. Fran the 
auto correlations, this energy is transferred into the lower frequencies. The 
transverse component auto correlations indicate little change in energy dis- 
tribution as the turbulence is contracted. The corresponding changes in the 
PSDs are shown in Figures 8 and 9. The change in the transverse structure of 
the turbulence as it is contracted can be found from the cross correlation, an 
example of which is shown in Figure 10. From the cross correlations, estimates 
of the transverse spatial auto correlations can be made by taking the values 
of the correlation at zero time delay and plotting them against the spatial 
separation of the probes. Examples of these spatial- auto correlations are 
shown in Figure 11 as the contraction ratio increases, and the decreasing 
transverse extent of the correlation area is apparent. From these plots, esti- 
mates of the integral length scales were made and examined in light of the 
analytical models in Section 6. The change in the mean square values of the 
component velocities is shown in Figure 12 and the initial suppression of the 
streemwise velocity and the amplification of the transverse velocity that is 
observed by other investigators (31,32) is also seen here. 

Turbulence fields are not the only type of distortion to be encountered in 
static engine operation. Steady and quasi-steady velocity distortions can also 
be present and the effect of a contraction on this type of distortion was also 
examined in the suction rig facility. The distortion generating bodies that 
were considered wore a 3.81 cm square rod, a synmetric NACA Series 6 airfoil, 
and a syitnetric NACA 0012 airfoil was used to generate a tip vortex. The axial 
measuronent locations were the same as for the turbulence measurements. 

The evolution of a steady axial velocity deficit as it is convected through 
the contraction is illustrated by the wakes of the airfoil and the rod shown 
In Figure 13. The airfoil wake is very rapidly and totally suppressed, while 
that of the bluff body persists scmewhat longer. 1h both cases no trace of the 
wakes was found at the -20.32 an (in the suction duct) location, at which the 
highest contraction ratios occurred. 
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Figure 12 Variation of RMS Velocity Through the Flat Contraction 


The NACA 0012 airfoil was set at the highest angle of attack possible without 
stalling, in order to generate as strong a tip vortex as possible. This was 
done by observing the flow over it with a cotton thread ...The angle of attack 
was about 14*^. The vortex was positioned approximately on the rig axis. Both 
streamwise and a transverse velocity component were measured in the vortex and 
the changes in these two aspects of the velocity field as the flow accelerates 
are shown ip Figures 14 and 15. The transverse, or azimuthal, velocity field 
is substantially unchanged by the flow contraction while the streamwise velo- 
city deficit associated with the vortex increases and broadens. This behavior 
of the streamwise deficit is completely opposite to that observed for the two 
dimensional wake and indicates the need for vortex suppression as near to the 
origin as possible. 
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5.1.? The Effect of Screens on Turbulence and Steady Distortions 

As ia the contraction experiment, the turbulence part of this test was con- 
ducted with the turbulence generator in place. Five screen types were examin- 
ed; three were made of honeycomb and two of perforated plate. The honeycomb 
used had a naainal cell diameter of .3175 cm and the thicknesses (cell lengths) 
were 0.953 cm, 1.27 cm and 2.54 cm giving thickness to cell dialieter ratios of 
3, 4 and 8 respectively. The corresponding honeycomb construction match al 
thicknesses were 0.008 cm, 0.005 cm, and 0.01 cm. The perforated plates that 
were tested had 41 percent and 66 percent open area. The hole diameter of each 
plate was .0366 cm and the plate thickness was .0126 cm. Each screen was 
mounted normal to the rig axis, 38.1 cm forward of the bellmouth, by means of 
a supporting frame and measurement locations were upstream and downstream of 
the screen at 30.48 cm and 48.26 oni (Figure 16). The duct velocities ranged 
from 30.49 m/s to 121.95 m/s. This variation allowed a range of screen 
throughflow velocities to be observed. A summary of test configurations and 
probe locations is shown in Figure 4. 

SCREEN 
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1 
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Figure 16 Axial Measurement Locations for Turbulence Screen Experiment 
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For the turbulence field defInUton the sane techniques were used as in the 
contraction experiment. Auto correlations PSDs, cross correlations and mean 
square values of velocity components were determined on either side of the 
screens. The effect of the various screens on the auto correlations of each 
velocity component is illustrated in Figures 17, 18, 19 and 20 with the cor- 
responding PSOs shown in Figures 21, 22, 23 and 24. The presence of small 
scale turbulence downstrean of the honeycombs is imuediately apparent. While 
sone of this may be residual frcm the upstream field, the majority is due to 
turbulence self-generated by the screen. This phenomenon is particularly ap- 
parent in the transverse tui-bulent velocity component which is. heavily sup- 
pressed by all honeycombs. The auto correlations (which are all normalized) 
downstream of the perforated plate, however, do not show a distinct bi-modal 
behavior which indicates that the self generated turbulence spans the whole 
spectrin rather than being confined to the higher wave ninbers. This is borne 
out by the PSOs downstream of the perforated plate, where the level Increases 
at all frequencies with increasing throughflow speed (Figure 25). The turbu- 
lent energy generation at the low frequencies is due to the unbounded nature 
of the flow and the finiteness of the screen. The high resistance perforated 
plates probably divert the incident flow over the edges of the plate where it 
separates and thereby produces turbulent kinetic energy at wave ninbers cor- 
responding to the plate dimensions. The PSOs of the velocity components down- 
strean of the honeycomb show the lew wave nmber energy to be constant, with 
the higher frequency levels increasing with throughflow velocity. The change 
in the transverse structure of the turbulence across the screen can be deter- 
mined to sone extent by the cross correlations. An exanple of these for honey- 
comb and perforated plate is shown in Figures 26 and 27. Estimates of the 
transverse integral length scale change are made fron these which are discuss- 
ed in Section 6. The streanwise interal length scales determined from the 
auto correlations are also exanined in Section 6. In view of the contamihation 
of the downstrean (post screen) turbulence field with self generated tur^- 
lence, it was necessary to extract the residual upstrean turbulence 
these results. This is discussed in Section 6 and in consequence the modifica- 
tion of the upstrean kinetic energy ccmponents by the screens is discussed 

there. 
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Figure 25 The Transverse Component PSDs Dcwnstream of the 66% Open Area 
Perforated Plate (Increasing Throughflcw Speed), 







For the steady distortion portion of this experiment, the generators used were 
a 3.81 cm square rod, a 5.08 cm x 7.62 cm rod and the afore mentioned NACA 
0012 airfoil for making a tip vortex. The rods produced sireamwise velocity 
deficits and the effect of the two screen types on them is shown in Figures 28 
and 29. In general, the higher the resistance of the screen, the greater is 
the suppression of the deficit. The screen having the highest resistance, the 
41 percent open area perforated plate produces a deficit inversion, changing 
the wake into a Jet. 


The tip vortex azimuthal vel ocity field was destroyed by all honeycomb tested, 
(Figure 30), however, the axial velocity distortion associated with the vortex 
increased across the honeycomb. It is interesting to observe that this distor- 
tion is almost cusp-like. The perforated plates were not effective as suppres- 
sors of the vortex, as can be seen in Figure 31. The core of the vortex be- 
cones larger due to the plates, but outside of this region, little change is 
seen in the azimuthal vel ocity f iel d. The axial velocity deficit change de- 
pended on the open area of perforated plate. The lower resistance plate in- 
creased it while the higher suppressed it. These contrary effects could be due 
to the previously mentioned deviation of the flow around the plate. 

5.1.3 The Pressure Drop Coefficients of the Screens 


Theoretical models of the pressure drop coefficients of perforated plate and 
honeycomb were presented in (29). Values of the pressure drop coefficients of 
the screens were used in the test program with a view to assessing these theo- 
retical models. Briefly the theoretical models are for a honeycomb, B1 as i us' . 
pipe flow coefficient formula. 
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While the results of the work of Baines and Pearson, discussed in 29, was sug- 
gested for perforated plate 


K 



1 

(1-s) 



2 


nr r • ■!= /•-! 


'Jr fs 


38 


oo 


96 CM MUN tVCOMB 


M% WtN AH E* PEWFOBM IPPi:^ 

VtlOClIYWf* 


VtlOCItVMrt 

»• 


o o o 


o 


o o 


o 


O o 

O Oh 

o 


DOVVNSTflGAM 

^ o o 




o o o 


o o 


o o 


’1 


untRCAM 

o O 


•? 


10 


—T — 
10 


transverse location cm 


rure 28 


The Charge in a Wake From a 
Through a Screen 


TRANSVERSE LOCATION CM 

3.81 cm square Rod on Convection 


,95 CM HONEVCOMO 


56% OPEN AREA PERFORATED PLATE 
velocity m« 


VELOCITY M/S 

*1 


8 


Q » 

□ 


a 3 


o 


□ 

o 


“ 1 

O o *o 


OOWNSTHEAM 

o o ^ 


.1 


□ 


o O 


* 


□ 


o 


i H 

O o 


o 


□ OOHHSTnIAM 


UPSTREAM 

O o o 


^■^^lrStRt AW 


•t 


-T" 

10 


4 0' 

TRANSVERSE LOCATION CM 




; 0 4 

transverse location cm 


~T~ 

10 


Figure 29 The Change in a Wake From a 5.08 cm x 
Through a Screen 


7.62 cm Rod on Convection 



AZIMUTHAL VClQpITYjj^ 


OOO 

ft 


velocity 
ft > 

o 


c o 

. 


<pL 


o O 


o 


o 


o 


STBtAMWIS e Vtl-OCI TY FIE ID 
v» locnv M« 

-O I- 


CM 






o o 


93 


s - 

c 

ft - 


o 


OOWNiTtUAM 


UKTREAM 

o O 

o 


-T r- “T T 

Xjj ft 0 » 

TRANSVERSE LOCATION CM 


10 


1 — r— T 

10 » 0 * 

TAANSVERSE LOCATION CM 


to 


Fijar^ 30 The Ctawe in the velocity Field of . Vortex on Conveotion Through 
1.2? cm Thick Honeycomb 


OD 


AZIMU THAL V ELOCITY Mt lO 
VtlOCITV M'S 


cP 


oq 


Si o' 


,U 


'1 O 


a 

a 


x> 


[J 0UlYI«Stm AM 
O iM*StlU AM 


o 


STREAMVTISE VELOCITY HELD 

VtlOCitvM/s 


O 

° o'* 


.q 


downstream 

O 


o 

O 


ftJ 


o 


>n 

CD 


UPSTREAM 

o O C) o 





TRANSVERSE LOCATION CM 



TRANSVERSE LOCATION CM 


Figure 31 The Chaise in the Velocity Field of a Vortex on Comection Through 
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The screen pressure drop coefficients were obtained by placing the screens 
over a nozzle and measuring the static pressure drop across them. This Infor- 
mation, together with the velocity at the screen (determined with a calibrated 
venturi) allowed the screen pressure drop coefficient to be calculated. A 
schematic diagram of the test arrangement Is shown In Figure 32. A gauze was 
mounted upstream in the supply pipe to produce a uniform velocity profile at 

the test screen. The pressure drop coefficient Is, by definition: 



PtU ” Ptd 3 

K = - 

or in terms of the measured quantities. 


1.115 CM NOZZLE 



Figure 32 Schematic Diagram of Screen Pressure Drop Coefficient Measurement 
Arrangement 
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These quantities were measured at several throughflow velocities and the re- 
sults are shown In Figure 33. The loss coefficient of each screen Is reason- 
ably constant with throughflow velocity. The higher solidity perforated plate 
(41 percent open area) has a higher resistance than the 66 percent open area 
perforated plate. The resistances of the 0.95 cm and 1.27 cm thick honeycomb 
screens are apparently reversed; however^ as was noted earlier In this sec- 
tion, the 1.27 cm honeycomb Is made of a thinner material. This raised the 
possibility that there was a significant base drag contribution to the honey- 
comb pressure drop. This possibility, however, was discounted- since the open 
area ratios of the screens were between 98 percent and 99 percent. On Inspec- 
tion of the a. 95 cm honeycomb under a microscope, It was seen that there was 
much more burring of this honeyccmb than the other two and the inconsistency 
In the pressure drop coefficients was ascribed to this. In Figure 34 a table 
is shown comparing measurements and theoretical predictions of the pressure 
drop coefficient. As can be seen for the 1.27cm and 2.54 cm honeycombs, there 
was a ratio of about 2:1 between measured and predicted values. The burred 
0.95 cm honeycomb had a corresponding ratio of 3.8:1.. Lt seems, then, that en- 
trance conditions to the honeycomb can be important in determining the pres- 
sure drop coefficient, K, and even with apparently smooth Inlet conditions 
there is still sufficient roughness to raise the value of K above the smooth 
pipe value. It should be recalled that in (29) the pressure drop coefficient 
of a 0.95 on diameter honeycomb was accurately predicted by Blaslus' formula 
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.3164 
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It is thus concluded that for any given honeyconb It Is possible for its pres 
sure drop coefficient to be up to four times the Blaslus' coefficient and a 
typical value of two is suggested for design purposes. 

2 ^ _ *6328 ^ 2 
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The measured values of the perforated plate pressure drop coefficients are 
consistent with the tneoretical values. Figure 33, and thus lend confidence to 

the model , 

5.2 The Transmission Loss and Directivity Effects of ICS Elartents 

In the acoustic part of the test progrem, the transmission loss of various 
screen materials in the presence of flow was measured. In addition, the effect 
on the radiation field of possible discontinuities in the ICS structure was 

examined. 

5.2.1 The Acoustic Transmission Loss of The Screens 

The screens tested in the aerodynamic part of the program were also used in 
the transmission-J-oss measurements. In the ICS configuration to be encountered 
on the test stand the angle of incidence of the acoustic wave to the screen 
material will be approximately normal . In any event, for the perforated plate 
this represents a worst case, and the honeyccmb transmission loss is very in- 
sensitive to incidence angle, (i.e. Transmission loss <ldB, for incidence 
angle <51°) 

The measuranents were made at P&WA's dual reverberation chamber facility. The 
source of noise was a Nor-Aire generator and flow through the duct was provid- 
ed by vacuun ^mps . A schematic diagrari of the test arrangement is shown in 
Figure 35. The test section lay between the reverberation chanbers and the 
duct was 2,54 cm in diameter to ensure only plane wave propagation over most 
of the frequency range of interest. The screen samples are shown in Figure 36 
and the method of mounting in the test section is shown in Figure 37. The no- 
minal throughflow velocities for which the screen transmission loss was mea- 
sured were 6.1 m/s. 9.15 m/s and 12.2 m/s. The noise generated by the Nor-Aire 
then propagates upstream through the screen. An array of three microphones was 
placed in each of the chambers. Figure 35. The sigials from each microphone 
were recorded simultaneous! y for a no-screen configuration, in addition to the 
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0.95 cm, 1.27 cm and 2,54 cm thick honeycombs and the 66 percent and 41 per- 
cent open area perforated plates. Subsequently, a one third octave band 
analysis was performed on the data from which estimates of the transmission 
loss were made. 


To determine the transmission loss of the various screen samples the noise re- 
duction between each chamber was first computed for each configuration. These 
noise reductions were obtained by subtracting the levels at each pair of mi- 
crophones for each one third octave band, i.e., 


Nn - SPLg - SPLj^ 


7 


Since there are three microphones in each chanber this operation produces nine 
estimates of the one third octave band noise reduction for each configuration. 
The transmission loss of each screen sample is then found by subtracting the 
no-screen (baseline) noise reduction from the noise reduction obtained for 
each screen sample. 

’’Sc ' '‘V - - '’"•R )sc - ( - SI’S ) B 8 

Equation 8 produces, therefore, eighty one estimates of each one third octave 
band transmission loss for each screen sarnple,_from which the average and 
standard deviation of the attenuation may be found. 

For the ideal situation of truly reverberant source and receiver chambers, the 
computations described above would be redundant since all microphones in any 
chamber at a given condition would produce the same measurement and one cal- 
culation would be sufficient to define a screen transmission loss. However, 
this is not the case in practice and the foming of many transmission loss 
estimates enhances the confidence in the result. 


The standard deviations of the one third octave band transmission losses for 
all screens and flows were, in general, of the order of the mean value. The 
majority of the values were then not significantly different from zero. 
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The average value of the 66 percent open area perforated plate transmission 
loss increases with throughflow speed. There Is however, little variation in 
that quantity with frequency, as can be seen in Figure 38. where the 9.14 m/s 
flow case is shown. At 12.19 m/s there is a considerable increase in transmis- 
sion loss, {Figure 39). The mean transmission loss of the 41 percent open area 
perforated plate however, decreases with increasing frequency at all speeds. 
With increasing speed there is small change in the average value, although at 
12.19 m/s. Figure 40, the scatter of the data is higher than for the lower 
speed conditions. The reason for the difference in behavior of the two perfor- 
ated plates in frequency is not understood. In addition, the higher transmis- 
sion loss measured for the higher open area plate at the higher one third oc- 
tave bands is intuitively questionable. 

The transmission loss of the 0.95 cm honeycomb is unchanged by increasing 
throughflow speed. The mean value is approximately zero for all one third oc- 
tave bands, except the highest where it increases to a maximum of about 1.5 
dB. The transmission loss shown in Figure 41 for the 9,14 m/s flow is typical. 
The 1.27 cr,. honeycomb has similar characteristics. The 2.54 cm honeycomb has a 
mean transmission loss that is more uniform in frequency than the thinner 
honeycomb, see Figure 42. In addition, the mean increased with flow speed from 
about 0 dB at 6.1 m/s to 2 dB at 12.19 m/s. In spite of the large standard 
deviation of transmission loss values encountered here, there is an indication 
that the thickest honeycomb (2.54 cm) may attenuate a normally incident wave 
at the highest flow speed {12.19 m/s). This suggests that the resistive com- 
ponent of the honeycomb impedance is no longer insignificant in this case. 

In surmary, then, there is a possibility that at the higher flow speed {12.19 
m/s) the perforated plate of 66 percent open area and the 2.54 cm honeycomb 
have non-zero transmission losses even though the standard deviations are 
large. In contrast, the thinner honeycombs {0.95 cm, 2.37 cm) present no 
significant obstacle to the incident sound field. The 41 percent open area 
perforated plate appears to attenuate the incident field in the lower one 
third octave bands. These results are compared with the theoretical predic- 
tions in Section 6.3 below. 
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Figure 40 
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Transmission Loss of 41% Open Area Perforated Plate at 12.19 m/s 
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Figure 41 Transmission Loss of .95 cm Honeycomb at 9.14 m/s 
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Figure 42 Transmission Loss of 2.54 cm Boneycomb at 9.14 m/s 

b.2.? The Effect on Acoustic Field Directivity of Honeycomb Corners and 
Structure 

The test progran was perfomed at Pratt Sf W hi tney Aircraft' s anechoic chamber, 
X-?07 stand. The test setup is shown in the schematic diagi'am in Figure 43 and 
in Figure 44. The chamber has a vol uae of 340 m and is lined on all sur- 
faces with specially constructed anechoic wedges to provide an anechoic 
envi niiiiiont for frequencies above 150 Hz. The floor, designed to be acousti- 
cally transparent, is constructed of tensioned wire mesh located above the 
f 1 oor wedges . 
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Figure 43 Schemdtic Diagram of Directivity Test Arrangement 



A metal frane fabricated from angle Iron was used to simulate one J 

be . Figure 45. Meeting at an angle of 135“. two .61m « .61m steets of 
63 cm tblc^ .95 on cell hone,conb were attached to the frane to fo™ 
angle of 135°. The corner fomed by the sheets of honeycomb was sealed with 

to pet sound leakage. The honeyc.b were renovable to allow fj test- 
‘n Of other corner structures, such as a 2.54 cm dianeter round rod and 
2 54 an square bar. Sound insulation material was added to all expose 
s'tructures (other than the test sections) to pre^nt spurious sound 
patterns. The sound source used was a 2.54 on dianeter »">vers,t. Driver Mo 
hor 1D65X located In line with the corner of the frane and the 2 

crophone as shown in ^ ^rectlvlty pattern. The driver was In- 

the driver to assure a smooth acoustic di recti v y v 

stalled on a rotating mount to all*, for changes of sound ’ 

50 10 “ incidence angles) with respect to the corner structure. The 

r;ver :ract1vated by a General Radio Model 1310 oscillator ^ ' 

tion with a McIntosh Model MC30 amplifier. The source was *iven noninal y 
2 kH 2 . 5 KHZ and 9 KHz. Capability to reorient the frane structure 
spect to the microphone array was also utilized. 
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(removed for BASELINE) 
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pigure 45 (f ometry 


of fjcneycomb Panel Mounting Frame 



Acoustic signals were detected by a pol ar array of .635 cm dianeter B&K- micro- 
phones (Model 4135) positioned at nonnal incidences to the sound source at a 
distance of 4.57 m. Microphone measurements were taken at 2° increments with 
respect to the sound source. Two microphone positioning methods discussed be- 
Icw were used to obtain the data; however, only Method (B) proved satisfactory. 

Method (A) 

Six microphones were located at angles of 110°, 115°, K0°, 125°, 

130° and 140° with respect to the test stand vertical axis. At each test 
condition the frame was repositioned angul arl y f rom 0° to 8°, in 2° in- 
crements, effectively providing measurements every 2°, from 102° to 
140°. This technique was curtailed during the test program when problans 
were encountered. The sound field within the chamber did not remain consistent 
for each position of the frame structure, resulting in a discontinuous direc- 
tivity pattern (i.e. repositioning of the frame produced a different sound 
field in frame fixed coordinates). 

Method (B) 


Twelve microphones were located every 2° from 110° to 132° relative to 
the test stand vertical axis. The angular range was necessarily smaller than 
Method (A); however, the resultant directivity pattern remained consistent. 
Test time was also reduced since fewer configuration changes were required. 

For each incidence angle tested a base line condition was obtained to define 
the directivity patterns of the driver and frame alone. To obtain the base 
line condition, sound levels were measured for each microphone at each of the 
nominal frequencies tested; 2000 H 2 , 5000 Hz, and 9000 Hz. Several repeats of 
of the base line conditions were made to assure measurement consistency. The 
measurements were repeated with the various test items attached; i.e., honey- 
comb, 2.54 cm s(^are bar and 2.54 cm round rod. For each of the microphones, 
the results were then subtracted from the base line data at the corresponding 
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mclOence ana fre<f,er.cy to give the airoctivHy wdiflcatton doe to the 

oresence of the body. The honeyconb corner configuration oas repeated at each 
Incidence angle with an alinilnun tape seal to determine whether sound leakage 
occurred with the screen cloth tape '.eal. The results Indicate no noticeable 
differences in directivit y patt erns. 

The directivity patterns obtained for the corner between honeycomb panels have 
been presented in (2g) and discussed here. An example of these patterns Is 
given In Figure 46 (The prediction shown Is that obtained from the theory e- 
scrlbed In (29)). In general more disruption of the far acoustic field was ob- 
served with Increasing frequency and Incidence angle. 

The two rods produced symetrlc directivity patterns and the modification to 
the free field pattern was also reater In amplitude at-the higher frequen- 
cies, as can be seen In Figure 47. The rod X-sectlonal shape does not affect 

1 the directivity pattern greatly. 


5.3 Conclusion 

The major observations from the test progran were 

0 The flow contraction amplifies the transverse turbulent velocity com- 
ponent and at low contraction ratios, it suppresses the stream-wise 
velocity conponent. At high contraction ratios the streanwlse turbu- 
lent velocity component Increases, although It did not attain the 

Starting value. 

0 Steady streanwise velocity deficits are suppressed rapidly by the 
flow contraction. The azimuthal velocity conponent of a streamwise 
aliped vortex was substantially unchanged by the flw contraction. 
The streamwise velocity distortion associated with the vortex in- 
creased in amplitude and extent. 
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Figure 46 


Directivity Bffeot of a aone»oomi Corner 





Figure 47 


acoustic Directivity Change Due to the Presence of a 2.54 cm 
eter Rod in the Radiation Field 
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0 The turbulence data gathered downstream of both honeycomb and perfor- 
ated plate was contaminated wUh self generated turbulence. The 
honeycomb screen suppressed the transverse turbulent velocity field. 
(The effect of the screens on other turbulence parameters is discuss- 
ed In the next section, where the residual upstream turbulence field 
Is separated from the screen self generated field). 

0 Steady streanwise velocity deficits are suppressed according to the 
resistance of the screen. In the resistance range studied here, the 
higher the resistance the greater the suppression. The azimuthal ve- 
locity field of the vortex was totally suppressed by the honeycmtbs 
and sotiewhat reduced and dispersed by the perforated plates. The 
stceamwise velocity distortion associated with the vortex increased 
on convection through the honeycomb and the 66 percent open area per- 
forated plate. 

0 The resistance of perforated plate is adequately predicted by theory. 
The resistance of the honeycomb is dependent on the entrance condi- 
tions and this precludes accurate estimates of the resistance for an 
arbitrary honeyconb. However an estimate may be made, by doubling the 
value obtained from B Iasi us' expression (Equation 6). 

0 The transmission loss of the thinnest honeycoabs (.95 cm and 1.27 cm) 
is negligible and invariant with throughflw speed. The thickest 
honeycomb (2.54 cm) has a transmission loss that increased with 
throughflow speed to 2dB at 12.19 m/s. 

0 The transmission loss of the 66 % open area perforated plate is negli- 
gible at the lower speeds (9.14 m/s and 6.1 m/s) but increases to 2dB 
at 12.19 n^s. This attenuation is constant in 1/3 octave band. The 
41% open area perforated plate transmission loss decreased with fre- 
quency and changes little with throughflow speed. The transmission 
loss is siyificant at the Iwer 1/3 octave bands. 
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„ The presence of hone/:«nb panel corners and structure can di sturb 
radiation field of a simple discrete source by slpnlflcant anoun s. 
Tim hl^er the frecmency the hl^er ttm disruption of the radiation 

field. 

The data gathered during the test progrm was used to assess and, 

s!^,. modify the theoretical models presented In (29). The results of this 

comparison are discussed In the next section. 
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6.0 THE DEVELOPMENT OF THE ICS DESICJ^ SYSTEM- {TASK G) 

In this section the test results of Section 5 are used to assess and/or modify 
the analytical models described in (29), Then the design criteria for an ICS 
are examined and the criteria for ICS location and material selection are es- 
tablished. The process is presented as a step by step procedure and exercised 
for the JT9D and JT15D engines. 

6.1 The Transfer Functions for the Effects of Contraction and Screening on 
Turbulence 

Measured changes in a turbulence field due to convection through 1) a flow 
contraction and 2) various screens, were described in Section 5. Theoretical 
models of these phenomena were identified in the Interim Phase II Report (29). 
The aim of this section is to compare the theoretical and experimental results 
and -thereby assess the models. If necessary the theoretical models will be 
modified. 

6.1.1 Contraction 

The comparison of measured with predicted characteristics as the turbulence is 
converted through a contraction begins at the coarsest level,. the kinetic en- 
ergy, and proceeds through increasingly detailed properties to the power spec- 
tral densities. 

The measured kinetic energyLratio across the contractions, as a function of 
contracti-on-ratio, is shown in Figure 48 together with the predictions using 
the simple models. The measured -vari at i on is approximately linear in contrac- 
tion ratio and its gradient is approximately half that of the theoretical var- 
iation. This discrepancy is in part due to viscous dissipation not accounted 
for in the model. A correction for viscous decay was suggested by Ribner and 
Tucker (33) and applied by Tucker and Reynolds (34) to their data with some 
success (i.e. their corrected data was in closer agreement with Ribner and 
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Tucker *^5 Theory). This correction is based on a linear decay of turbulent kin- 
etic energy with distance in a unifonn mean velocity field, and can be written 
as 


where in Equation 9 refers to the position of the virtual source of the 
turbulence. This correction was applied to the present data and the corrected 
values are shown in Figure 49. Also shown in this figure are data gathered by 
Uberoi (31) and Hussein and Ramjea (32). These data are also viscous decay 
corrected and as can be seen, the data are in good- agreement with each other. 
Various contraction types were studied by these investigators; Uberoi examined 
contractions of square cross section, Hussein and Remjee considered circular 
cross sections. The contraction rates were also varied. The agreement is quite 
close and thus the decay process through the contraction is not strongly de- 
pendent on the details of the contraction. 


On comparing the viscous decay corrected data to the predictions of the 
Satchel or-Proudnan, Ri bner-Tucker theory, much closer agreenent is observed 
(Figure 50). The simple theory predictions are not as close as the Batchelor- 
Proudnan, Ri bner-Tucker theory. The corrected data is still sonewhat low at 
the-higher contraction ratios, but this is attributable to error in estimating 
the location of the virtual turbulent source, i.e. X 


The next higher level of detail of the turbulence field is supplied by the RMS 
values of the component turbulent velocities. Two vel oci ty components were 
measured in the test program and the resulting RMS velocity ratios through the 
contraction are shown in Figure 51. These data are not viscous decay corrected 
in this figure. The predictions using the simple theory of Prandtl discussed 
in (29) overestimate both the increase in the transverse velocity component 
and the suppression of the streamwise velocity ccmponent. The theory of 
Satchel or- Proudm an, Ribner-Tucker produces predictions that are closer to the 
data but still err in the same way as the simple theory predictions. A reason 
for the departure from the linear theory (in addition to the viscous decay 
discussed earlier) was foretold in the Interim Phase II Report (29). There it 
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Figure 48 Kinetic Energy Ratio Through the 
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Figure 50 The Kinetic Energy Ratio Through a Contraction Shoeing the Effect 
of Correcting for Viscous Decay 



Figure 51 The Effect of a Contraction on the Turbulent RMS Velocity Ratios 



was stated that for this configuration the turbulence inertial effects are not 
negligible and that consequently there would be a tendency to return to iso- 
tropy in the contraction (or perhaps more accurately, a tendency to equiparti- 
tion of the kinetic energy between the three velocity components). This ten- 
dency is seen most apparently in the streamwise component (Figure 51) as it 
reaches a minimun and subsequently increases. This type of deviation from the 
linear theory is not explainable by viscous decay. The tendency to equi parti- 
tion is shown more clearly in Figure 52 where the ratio of the RMS values of 
the transverse to streamwise velocity component is plotted as a function of 
contraction ratio. At the higher contraction ratios a maximum of this RMS ra- 
tio is reached and subsequently a tendency to equipartition of kinetic energy 
is seen. The predicted values are also shown on this figure. The anisotropy 
predicted by the linear theory is observed in the data only at low contraction 
ratios ( <2). In this regime for the present configuration this probably cor- 
responds to the region where viscous decay and turbulence inertial effects are 

negligible. The data of Uberoi and Hussein and Ramjee will now be compared to 
the current data. 

Figure 53 compares the data of Uberoi and, Hussein and Ranjee. The first thing 
to note is that the close agreenent in kinetic energy ratio data of the three 
experiments, previously noted, is repeated when the transverse RMS velocity 
ratios are compared (i.e. The ratios obtained by dividing the RMS value at a 
contraction of 1 by the RMS value at a contraction ratio of unity). The 
transverse RMS velocity ratio data of Uberoi agrees very closely with the pre- 
sent data. The data of Hussein and Ramjee is somewhat lower than Uberoi 's and 
the current data, but those authors note that the position at which their ref- 
erence RMS values were determined was not the same as Uberoi 's (who used the 
beginning of the contraction). !!ussein and Ramjee assert that their transverse 
RMS velocity data agrees closer to that of Uberoi if the same reference posi- 
tion is used. Thus the variation of the RMS ratio of the turbulent transvecs^ 
component is to a large extent independent of the contraction rates and shapes. 
The viscous decay correction for the three experiments is similar, so that the 
same statement can be made of the corrected transverse RMS ratio data. 
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COMt RACTION RAI10 ' i . 

Figure 52 Ratio of Transverse To Stream/ise RMS Values Through The Contraction 
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Figure 53 Comparison of Present RMS Ratios With The Data of Uberoi, And 
Hussein And Ramiee 


In contrast, the RMS ratios (i.e. Ratio between the RMS value at contraction 
ratio and at a unity contraction ratio) of the streamwise turbulent 
velocity component for the various contractions are quite different although 
the general trends are similar. This appears to be in conflict with the pre- 
vious observation that the total kinetic energy ratio and the transverse 
velocity ratio variations are approximately the same, for the experiments of 
Uberoi, Hussein and Ramjee and the present data. For if the kinetic energy is 
the sanje and the transverse components are the same, one would expect the 
streamwise conponent variations to be the same for all three experiments. How- 
ever, if the streanwise turbulent velocity contribution to the kinetic energy 
is small compared to that of the transverse components then variations in the 
streamwise velocity component level will not affect the magnitude of the kine- 
tic energy significantly. This is the case for the contractions considered 
above, except at the lowest contraction ratios. 

In view of the dominance of the energy in the transverse turbulent velocity 
components, it is important that the behavior of this component be properly 
modelled. Thus the viscous decay correction is applied to these components. A 
comparison of the corrected transverse RMS levels with theory is shown in 
Figure 54. There is close agreement with predictions using the Satchel or- 
Proudnan (35). Ribner-Tucker theory. The simple theory is not as good. The 
Satchel or-Proudnan, Ribner-Tucker theory when corrected for viscous decay then 
provided the best estimate for the transverse RMS ratio through a contraction. 


Neither this theory nor the simple theory are close for the streamwise RMS. 
ratio predictions although at the low contraction ratios (<6) before energy 
transfer is apparent, the Satchel or-Proudnan. Ribner Tucker theory provides a 
reasonable estimate. The model for this component will be selected based on 
this and other considerations to be discussed in Section 6.4. 
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The in length tele reties (I.e. The ratio between t e '» 9 t c le 

at a rontractlon ratio of «i and the value at a contraction ratio of unity) 
through the contraction are shown In Figures 55 and 56. Four Integral long h 
scales were computed from measurement {as described In Section 5). an axia 
and a transverse scale of the streamwise turbulent velocity component (Lii. 

C) and an axial and a transverse scale of a transverse turbulent ve ocity 
onf n L ^ Considering the data first, it is interesting t 

component (Ljj, jj). , of energy observed between the turbu- 

note that the tendency to equlpartition or eneigx v.. 

Teli! v!loc1ty components Is not reflected by a tendency to Isotropy In t e 
integral length scales. The axial scales Increase and the transverse scale 
decrease almost monotonically. Both axial scales and 1° , 

of the streamwise turbulent velocity component, ’* 

predicted by the simple theory. The scale ratios of the transverse component. 
(I.e. ^21B ). however, are not well predicted at the higher 

^21A ^23A ^13B , ^23B 1 

contraction ratios. The measured component scale ratios, (i.e. _ ' j— 


are closer to the predictions of the viscous-decay corrected Ribner-Tucker 
theory (33). also shown In Figures 66 and 56. The length scale ratios were 
determined using this theory by U. Bant (36). He found that for the transverse 

velocity components 
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Figure 55 The Axial Length Scale Ratios Through A Contraction 
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Figure 56 The Transverse length Scale Ratios Through A Contraction 



The predictions in Figures 55 and 56 used the expressions 11 and 12 with the 
viscous corrected mean square velocity ratio discussed earlier in this sec- 
tion. In consequence the ICS design system will use the simple models for the 
streanwise component length scales andGanz's formulation (with the viscous 
corrected mean square velocity ratio) for the transverse velocity component 
length scale ratios. 

Two types of PSDs of the turbulent velocity field were determined by measure- 
ment, (see Section 5), both in the frequency domain (which corresponds to the 
streamwise wave number under the frozen turbulence assumption) of the trans- 
verse and streamwise velocity components. The pre and post contraction PSDs 
are shown in Figure 57 for the transverse velocity component. There is very 
little change in this PSD over the contraction range examined. It should be 
noted that in this figure and Figure 58 only the spectral shapes are being 
compared (each PSD has been normalized by it's maximun value). This result, 
i.e., the constancy of the transverse component PSD, is inconsistent with the 
simple model since from the measurement 


where F denotes a normalized PSD and noting that 

(0 = kj^g 15 

This is Taylor's frozen turbulence assumption, and the PSDs in this equation 
and in the subsequent section, show the frequency domain may be mapped into 
either the upstream or downstream axial wave number space using these trans- 
formations. 


Then Equation 14 becomes 
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Normalized Streamise Velocity Component PSPs At Various 


Contraction Ratios 
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This Is In conflict with the simple model which would predict that 
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The-approxiraatton of Equation 16 is. houever predicted by the Ribner-Tucker 
theory (33). In their report the comparison of the PSD, Fj (k^) before an 
after tlie contraction shows that the PSOs are very similar in shape when plot- 
ted as a function of the upstream streamwise wave number. 


The streamwise velocity component PSOs are shown in Figure 58. Each PSD is 
denoted by the contraction ratio at which it occurs. Those marked with 1 are 
pre contraction PSOs. If this figure is viewed together with Figure 51 it is 
iimiediately apparent that as the streamwise component RMS value becomes lower 
its spectrun shape diverges from the pre contraction PSD. At a contraction 
ratio of 5.65, a minimun in the RMS value is reached corresponding to a PSD 
shape that is most different from the pre contraction PSD. Subsequently as 
energy is transferred into the streamwise component from the transverse com- 
ponent. the streamwise PSD becones more similar to the pre contraction PSD. In 
the pre energy transfer region (i.e. contraction ratios <5.65). if the PSOs 
are plotted in terns of the local wave mmber, it is clear that they are very 
similar in shape. An example is shown in Figure 59 where the PSD at a contrac- 
tion ratio of 2.29 is cmpared with the pre-contraction PSD. In short, for 
contraction ratios <5.65. 
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This is in accordance with the simple model. The characteristics of this PSD 
as predicted by the theory of Ribner and Tucker (33) and in the extension o 
U. Ganz (36) are not observed In the streamwise velocity canponent data. 
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From Equation 19 and the corresponding mean square velocity and Integral 
length scale ratios, the relationship between the un-nonnallzed PSDs may be 
determined as 



WAVENUMBER - m '' 
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For the other streamwise velocity component PSDs, it Is reasonable, in view of 
the good agreement of the predicted length scale ratios with measuranent , to 
assure that these PSDs will be related in a similar manner to Equation 20. In 

general then for the streamwise velocity component in the initial contraction 
region: 



At the higher contraction conditions the streamwise component spectra change 
radically. Tnis is due to the transfer of-- energy from the transverse velocity 
corponents discussed earlier (recall that the energy transferred from the 
transverse velocity component is small initially compared with the total 
transverse velocity kinetic energy, but is large, compared to the -existing 
streamwise kinetic energy. Therefore while the effect of this transfer. on the 
transverse spectra is small, the effect on the streamwise velocity component 
is large). Specifically the streamwise PSD receives energy in the low fre- 
quency regime and this is expected since the mechanism of transfer is iner- 
tial. Neither the simple model nor the linear theor^ of Ribner and Tucker pre- 
dict the effect of energy transfer on the streamwise conponent of the turbu- 
lent velocity field. However, there is an alleviating condition that will be 
discussed below. 

The comparisons of measured and predicted PSDs together with the scale and RMS 
velocity ratio conparisons point to the conclusion that the transverse charac- 
teristics are most accurately portrayed by the viscous decay corrected Ribner- 
Tucker, Batchelor-Proudnan theory. Thus the 6 transverse velocity PSDs will be 
modelled in the ICS desiyi system by the viscous decay corrected formulations 
of Ribner and Tucker and Ganz (33,36) who extended the work of Ribner and 
Tucker to complete the PSD matrix. The streamwise velocity scales and PSDs in 
the pre-energy transfer region, described above, are accurately modelled by 
the simple theory. Other characteristics, namely, the RMS vel ocl ty rati os and 
PSDs in the energy transfer region are not described well by either theory and 
it is app*-opriate to discuss at this time the alleviating condition mentioned 
above that lessens -the importance of accurately predicting streon wise charac- 
teri sties . 
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The coinponent of the turtxilent velocity field that is probably most important 
as regards noise generation is the upwash in the tip region of the fan rotor. 
As shown in (29) this is related to the transverse and streamwise turbulent 
velocity components by means of the relative blade inflow angle, “Iq , 

(i.<) (•“) <■'<'»«"' Vq *■ ^''2 

where F denotes a PSD and the subscripts 1 and 2 refer to the streamwise and 
transverse velocity components, respectively. 


Integrating this equation over frequency yields, 

— 5 ^ 2 ’ 2 2 

cos'’ Vq + U2 sin ^ 
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Now if it is assimed that the turbulence is initially isotropic then 

and the ratio of the mean square values of the upwash com- 
ponent across a contraction may be written as 





sin > 
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where B denotes the post contraction location and A, the pre- contraction, or 
upstream location. In the tip region of a JT90, a typical value of is 
30°. If this fact is used together with the measured data at a contraction 
ratio of 21.3 then an estimate of the mean square ratio of the upwash velocity 
component across the contraction can be made, and is 



.2Bin + 2.48 
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The upwash velcx:1ty ratio thus appears to be dominated by the transverse tur- 
bulent velocity component and hence the prediction of the turbulent upwash 
velocity through the contraction Is probably not strongly affected by errors 
in the streamwise component characteristics provided the energy transfer Is 
not too far advanced. Consequently for the purposes of the ICS design system 
and in view of the tendency of the streamwise catiponent energy to recover at 
the high contraction ratios encountered on the test stand (see Figure 53, 
where at a contraction ratio of 21 the streamwise RMS ratio c? .6), It Is sug- 
gested that the streamwise component mean square velocity ratio be approxi- 
mated by .5 In the higher contraction regime 6, say). With regard to 
the post contraction streamwise PSD, In this higher contraction region, there 
Is some similarity between this PSD and the pre contraction PSD in the up- 
stream wave number, It Is thus suggested that the PSD relationships be 
at the higher contraction ratios, for purposes for ICS design. 



In suimiary the contraction models to be used In the ICS design system are 
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6.1.2 Screening 


The experimental determination of the effect of a screen on a turbulence field 
conuectlng through It, Is complicated by the turbulence generated by the 
screen Itself. This screen- generated turbulence contaminates any downstream 
turbulence measuranents as noted In Section 6. Furthermore since the energy 
lost from the mean flow Is converted In part Into turbulent kinetic energy by 
the screen, the greater the mean flow energy loss the greater will be the tur- 
bulent kinetic energy gain. Thus the screen generated turbulence Intensity 
level can be expected to Increase with Increase In screen resistance and 
throughflow speed. 

The cent an i nation of the downs trean turbulence field was particularly evident 
in the perforated plate configurations. This contamination is illustrated in 
Figure 60 where pre and post screen transverse velocity component PSDs are 
shown for various throughflow velocities for the 66% open area perforated 
plate. As the duct speed increases (this is the velocity in the suction rig 
duct Figure 5), the throughflow speed increases and the self generated turbu- 
lence increases. It should be noted that the PSD level increases at all fre- 
quencies. The high frequency increase is due to turbulence generated at the 
perforated plate holes while the increase at lower frequencies is probably due 
to the large scale turbulence generated by the flow arouncLlhe edge of the 
screen. This latter source of self generated turbulence will not he present in 
an ICS. The small scale hole-associ ated turbulence, however, will be present 
but will not contribute to the fan tones provided the perforated plate hole 

si ze is small . 

In order to extract the characteristics of the residual upstream turbulence 
field, uncontaminated by screen generated turbulence. It was felt that since 
the downstrcOT Intensities and kinetic energy Increased as the throughflow 
increased. If the data were extrapolated back to small throughflow velocities, 
a reasonable estimate of the residual turbulence characteristics would be ob- 
tainable. In Figure 61 the kinetic energy ratios across the perforated plates 
are shown as firctlons of duct velocity. On extrapolation back to zero duct 
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I (this corresoonds to a low throughflow velocity), values of the ktne 
speed (thi turbulence fields are obtained for 

energy ratio due prim y extrapolation Is 

the two perforated Te iron throughfl» velocity Is not as 

perfonned using „,th simple model predictions and 

accurately known. These da » . ca co thp aareement is 

Townsend's (37) data (obtained using gauzes) In Figure 62. The agreement 

good. 

j a. .. 4 -sms.i-o tho Tstlos of thc Hiean square values 
Figure 63. As can be s , .(.otroDV for each perforated plate (it 
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different from the theoretical pre ic ion. extent of the 

es^pir:" retum to 

unequal energy distribu ion u 

Isotropy provides the lower bound. The choice o tUuatlon In 

velocity ratio will be discussed in the context of the test 

Section 6.4, below. 

ratios were determine for each perforated plate. They are ll- 

L u When. for Instance. Is the Integral length scale of 
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the u, velocity component in the Xj direction. 
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with either theory. This Is not 

The measure.) scale ratios ^ ^ the n>e.o s,e«-e veloclt, ratio 

enexpected In view of the to the retornln, of the turhulerce 

to isotropy since this ^tlon experiment, while the energy dls- 

valoe. Thus, as observed In the there Is not a correspond- 

trlbutlon tends to ftHttaln their Isotropic values. The scale 

lug tendency for the length sea to the diffusion of the flow 

data is. however, to sane extent gues ^,^tlc en- 

the perforated plate. This dlffus^ Conseguently the scale model will be 
ergy ratio across the tHles discussed below In sub section 

chosen as a result of other velocity exponent downstream of 

6.1.3. The shape of the PSD o frequencies 

the perforated plates Is s m -nncles are contaminated with self- 

(Figures 66 and 66). The higher requ accorda,.oe with 

generated turbulence. This ' jf,geres 67 and 68) tend to 

the simple theory. The streamwlse gre-sUen PSO's as the plate 

become similar In shape to t e cor ^ 5 ,„p,e 

throughflow velocity decreases. Thus th« terbulence Is unchanged 

theort that, as a whole. -;t:rrtprr::;;of the tendency to return to 
by the perforated p a e. ..re velocity components. 

Uotropy observed in the mean square velo 

.incltv ratio's indicate a tendency to 
,n sumnary. both PSD's and 7Jth''scales do not yield clear Information 
return to Isotropy. The Integr^ plate Is predicted fairly well 

The kinetic energy change ac turbulence PSD shapes downstream 

by the simple theory as je th ,„,P„,e„ce w 

the perforated plate. The effect 
now be considered. 

f .11 honeycombs examined was also contmlna- 
The turbulence field downstream o the PSDs 

ted by self-genorated particularly noticeable In the auto cor- 

::ratir:-:rT,re:::’v:-ny conponent. Figure 6,. Here the Increasing 
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level of high frequency turbulence Is seen with Increasing through flow velo- 
city. This Increasing level of high frequency energy is manifested in the 
autocorrelation as a progressively sharper spike. The upstrean autocorrela- 
tions are also shown In the figure, and the higher frequency energy content- of 
the downstream turbulence field Is self-evident. The extractiwi of the resid- 
ual turbulence levels and Integral length scales Is carried out by using the 
auto correlations, (see Figure 70). (This procedure Is used rather than the 
extrapolation method described for the perforated plates because It produced 
more data points and also because the self generated turbulence was not clear- 
ly defined In the perforated plate auto- correlation.) Thus the turbulent kine- 
tic energy ratio across the honeycomb could be found for the turbulence con- 
vected through the honeycomb, uncontaminated by self-generated turbulence. (It 
is unlikely In view of the small scale of the self generated turbulence that 
It will contribute to fan tone noise, however,. there. could. be a contribution 
to fan broadband noise). 

The ratios determined from the residual convected kinetic energy are shown In 
Figure 71, They are plotted as a function of the pressure loss coefficient. 
This paremeter Is chosen rather than the flow angle ratio («) which Is effec- 
tively zero for each honeycomb (see Section 5). The theoretical values of the 
kinetic energy ratio are also shown as a function of pr essure loss coefficien t 
for O' = 0. 

^‘ 1 . - i / 1 \ ^ 

KE^ “ 3 \T+k/ 40 

The agreement with the predi ction of the simple theory is good. 

The mean square component velocity ratios (with the self generated contribu- 
tion removed) are shown in Figure 72. The streamwise velocity component levels 
are higher than the transverse component levels. The equi partition of energy 
observed downstream of the perforated plate Is not therefore, repeated In the 
case of the honeycomb, as might be expected, since honeycomb greatly attenu- 
ates the transverse velocity component. On comparison with predicted levels, 
however. It Is noted that the measured transverse component levels are higher 
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than the predicted values (of approximately zero) while the streamwise compon 
ent values are for the most part, lower. A slight tendency to Isotropy or 

equipartition of the turbulent kinetic energy: rel ative to the theoretical pre- 
diction Is thus indicated. 

Again then, as in the case of the perforated pi ate, the data lie between val- 
ues predicted assuming a return to isotropy and those predicted assuming no 
return. The situation likely to be encountered downstream of an actual ICS 
will be discussed in the following section. 
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f’igur#' 72 Cof7ip,»r ison ,.f Nr.isurt-J an,f Pirdirted MfM/> Velocitu Ratios 

Across Hoiit'Qcomb 


The Integral length scale ratios are shown in Figure 73. The transverse scale 
ratios tend to be close to the predicted value of unity, while the axial 
length scale ratios are on the whole somewhat higher. The 2.54 cm thick honey- 
comb (K “ .79) yielded erratic values and this was due in part to the low 
level of the residual turbulence compared to the honeycomb generated level. 

The length scale ratios support the non-isotropy in the turbulence field ob- 
served in the energy distribution. The length scale ratios may be roughly 
approximated by unity and this value will be used in the ICS design system. 

The change in shape of the PSDs of the two velocity components is shown in 
Figures 74 and 75. The 2.54 cm thick honeycomb configuration is shown here, 
but the comparisons are typical for the other honeycombs. Each PSD is normal- 
ized by its maximum value thus allowing a comparison of the shapes to be made. 
The comparisons of the streamwise velocity component PSDs, Figure 74^ indicate 
that the honeycomb generated turbulence increases with increasing throughflow 
speed. This corresponds to the deductions made previously from the autocorre- 
lations. In addition there is a tendency for the downstream PSD shape to ap- 
proach the upstream PSD shape as the throughflow velocity decreases (i.e. at 
the 30.48 m/s duct velocity condition, corresponding to the lowest throughflow 
speed, the downstream streamwise velocity PSD is most similar to the upstream 
PSD shape). The transverse component PSD comparisons show a similar trend, 
however, because of the extreme suppression of the component, the downstream 
PSDs are dominated by the self- generated turbulence even at the lowest through- 
flow velocity. The current measurements then provide little information on the 
distribution of the residual energy in the transverse velocity component down- 
stream of the honeycomb. The integral length scale ratios discussed above are 
roughly unity, which suggests that the residual energy distribution of the 
transverse velocity component could be similar to the prescreen distribution. 
These data then do provide some support for the simple model PSD relationship 


where? represents a PSD normalized with respect to its maximun value. This 
model will then be used in the ICS design syston. 






POST SCREEN 



In sunn ary, the turbulence kinetic energy ratio across both types of screen 
(perforated plate and honeycomb) with self generated turbulence corrected out 
Is reasonably well predicted by the simple theory of Taylor and Batchelor. The 
corresponding values of the streamwlse and transverse components and mean 
square velocity ratios are not well predicted for the perforated plate In view 
of an apparent return to Isotropy. The theoretical predictions are closer for 
these ratios across honeycomb, probably because the honeycomb produces a much 
more uneven energy distribution between streamwlse and transverse velocity 
components. The measured Integral length scale ratios across the screens are 
scattered around unity, while the downstream PSOs are contaminated to a 
greater or lesser extent by self- generated turbulence. There Is sane evidence 
to suggest that the residual turbulent energy Is distributed In wave mmber In 
a similar way to the pre-screen turbulence. 

The models for the effect of perforated plate and honeycomb on turbulence to 
be used in the ICS design system are; 
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6.2 The Transfer Functions for the Effects of Contraction and Screening on 
Steady Distortions. 

As described previously (Section 5) three types of wake generators were tested 
to study the effects of contraction and screening on steady distortions. Two 
of these produced streamwise velocity deficits while the third, an airfoil, 
produced a vortex. For the contraction part of the program a 3.81 cm square 
rod and a symmetric airfoil, 1 cm thick, were used to produce the axial 
velocity deficits. For the screening part., a 3.81 cm square rod and a 5.08 cm 
X 7.62 cm rod were used to produce the axial velocity deficit. The data 
gathered in these experiments were used as described below, to validate or 
modify the theoretical transfer functions. 

6.2.1 Contraction 


The transfer function for the effect of contraction on a steady axial deficit 
is first considered. From the measurements using the 3.82 cm square rod, 
values of the wake deficit ratios and wake width ratios as a function of con- 
traction ratio were obtained. The tabulation of these ratios is shown in 
Figure 76. The comparison with theory is shown in Figure 77. Only two values 
of the measured wake width ratio are shown in Figure 77, since at the other 
contraction ratios the velocity deficit ratio was zero and the corresponding 
wake width ratio was indeterminate. 


The theoretical model for the contraction of a viscoiis-jwake was reported in 
the interim Phase II report (29) of this contract. 
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Figure 76 Wake Ratios Through a Contraction (3.81 cm Square Rod) 



CONTRACTION HATIO ' (, 

Figure 77 Comparison of Measurement and Theoretical Prediction (3.81 cm 
Square RodJ 





















The value of the virtual source distance X(r^j) was, for a case with the con- 
traction ratio equal to unity determined from the meas.ured data and was found 
to be .29 cm. The model predictions In Figure 77 use this value. 

The comparison of measured with predicted values shows that the trend of wake 
velocity deficit decay is duplicated by the model but the absolute velocity 
deficit ratios are over predicted. At low contraction ratios, not unexpectedly, 
the agreaaent Is close. The wake width ratio is not well predicted by the 
model. The two data points available indicate an Increase in wake width ratio 
with contraction ratio. However, it Is expected that at high contraction ra- 
tios the slope will be negative. A mechanism not accounted for. by the theory 
is thus indicated. 

The discrepancies between model predicted ratios and measured values could be 
due to an interaction between the contraction and viscous effects unaccounted 
for by the model. A possible explanation which would explain the data is that 
at low contraction ratios the viscous forces in the wake are more dominant 
than the inertia forces thus causing the wake width to Increase. The velocity 
deficit correspondingly decreased more rapidly than would be predicted by a 
theory not accounting for this interaction. Subsequently as the contraction 
r-atio increases, the inertial forces become increasingly dominant causing the 
wake width to contract. 

To account for this contract! on- viscous interaction the following procedure 
was-adopted; 

a) The distance between measuring stations D was adjusted to fit the measured 
velocity deficit ratios. With this approach the interaction effect mani- 
fests itself as an effective distance D which is a function of contraction 
ratio, In normalized form, using the wake generator thickness, a, 
this function was determined to be 

^-30.12 t, - 21.56 

o 1 ■ cn 


and is shown in Figure 78. 
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CONTRACTION RATIO ~ V, 


Figure 7B Normalized D-Bffective as a Function of Contraction Ratio 


The corresponding value of the virtual source distance X(r ), when normal- 
ized is ° 

X(ro) 


b) A multiplicative constant was applied to the wake width ratio expression 
and evaluated using the measured data interpolated at i ■ 1. This con- 
stant was determined to be .318. 

The adjusted transfer finctions thus become: 

Velocity AU(r) _ „ -1/4 ) 3/2 (30.12 a. - 21.56) n, / 

Deficit AU(rJ ~ h V^l + — > 


Tlie plots of these models are shown in Figure 79 together with the measured 
val ues . 

From the measurements obtained using the s>mnietr1c airfoil, shown in Figure 
80, it is apparent that the wake generated by this type of body recovers very 
rapidly. In fact at the 25.4 cm and -20.3 cm measuring stations at which data 
at contraction ratios ranging fran 1 to 14 approximately were gathered, no 
wake was detected. In the plots of Figure ao the predicted values of the ve- 
locity deficit and wake width ratios using Equations 52 and 53 are shown. Ex- 
pression 51 above was used to determine the virtual source distance. A lew 
value of the velocity deficit ratio is predicted using this model. 

While the adjusted expressions fit the measured data well there is no guaran- 
tee that these expressions are universally applicable. However, these measure- 
ments have indicated that a viscous wake in a contracting flow is subject to 
an accelerated recovery in excess of that expected by a linear combination of 
viscous and contract! on- effects. With this in mind and considering the model 
predictions of the velocity deficit ratio at high contraction ratio. Figures 
79 and 80, it seems that the velocity deficit as predicted by pure contrac- 
tion AU(r) -1 provides a valid upper bound. Similarly the 

AU (T^T " '^1 

prediction of the wake width provides a minimum estimate ot_ 

that quantity. 


The effect of contraction on a vortex will now be considered. 


A vortex is comprised of two regions, an outer region where viscous forces are 
negligible and an inner core region where they are not. When a vortex is con- 
vected through an axisymmetric contraction, consideration of the conservation 
of circulation leads to the conclusion that the outer velocity field is un- 
changed, see Appendix B. In the vicinity of the core the theoretical model 
discussed in the Interim Phase II report (29) predicts that the maximim trans- 
verse velocity increases according to 
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Figure 79 Comparison ot Adjusted Model Predictions With Measurement (3.81 
Square Rod) 



Pigure 80 Comparison of Adjusted Model Predictions With Measurement 
fSummetric Airfoil) 



and the core 


radius reduces according to 
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The transverse velocUies measured at the three axUI locations 48.26 on.. 
30.48 cm., and 20.32 cm. (referred to the bellmouth) through the 
s! port t^e invariance of the transverse velocity field In the out^ r Ion 
(Figure 13). The maximmn velocity (at the core boundary) changes slightly. 

in the light of these data It Is concluded that the viscous Inertial force 
egulllbril at the core boundary of the vortex Is not significantly changed y 

convection through an axl-sjmmetric flow contraction up to » 

of about 15, for the contraction lengths examined here. In fact the core - 

Plus of the vortex renalns relatively unchanged (Figure 81) as to 

tion and therefore the vortex transverse velocity field. Accordingly the ana- 

lytically-deri ved functions are modified 
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Figure 


81 The Variation of Vortex Core 


Radius With Contraction Ratio 


It w«s obsarved In Section 6. that there is an axial velocity distortion asso- 
ciated with the vortex. This distortion is roughly a deficit. There is ev - 
dence (Figure 15) that this axial velocity deficit is cusp shaped at the be- 
ginning of the contraction. Subsequently as it is convected into 
tion the deficit increases and becomes fuller and broader (Figure ). 
deficit does not then behave as a two dimensional wake. Figure 82 presents the 
variation of axial velocity deficit ratio with contraction ratio. It should^be 
noted that at the initial part of the contraction it is difficult to detemine 
the minimus of the axial velocity due to the cusp nature of the deficit. M 
the higher contraction condition measurement of the full recovered axial v^ 
locity component were Impeded by the presence of the duct. Ml deficit esti- 
mates in this Figure (82) are therefore probably low. Mso plotted in this 
figure is the theoretical axial velocity deficit ratio for an inviscid wake 
for reference. AS can be seen initially (at low contraction ratios) the mea- 
sured velocity deficit decreases (data follows approximately this curve), but 
as the contraction induced inertial forces Increase, the deficit on an abs^ 
lute basis tends to Increase. The broadening of the deficit in this phase is 
probably due to viscous diffusion. Thus viscosity disperses the axial velocity 
distortion more rapidly than the transverse (vortical) velocity distortion. 



Wpure ei v.irlation of Axial velooltg oefloltg Ratio in a vortex with 
Contraction Ratio 
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From a noise generation standpoint the transverse- velocity field of a vortex 
Is Invariant under the axisynmetric contraction of moderate length, while the 
axial velocity distortion associated with the vortex Is amplified on an abso- 
lute basis by the contraction. This Indicates that If a vontex exists In the 
contracting Inflow to a fan on a static test stand. It should be suppressed 
(e.g. by use of an ICS) as far away from the fan face as possible. 

In suTsnary the steady distortion contraction models to be used in the ICS - 

Design System are 


1) Axial V el ocity Distortion 
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6.2.2 Screening 


As In the previous discussion on contraction defects, the modification of an 
axial velocity distortion will first be examined. Five screens were used and 
as found in Section 5 the pressure drop coefficient, K, varied from .41 for 
the thin honeycomb to 5.9 for ti. 41% open area perforated plate (See Figure 
33). The generators of the axial velocity deficits were the 3.81 cm square rod 
and the 5.08 cm X 7.62 cm rod. For each rod and screen conbination, three duct 
conditions were set at 53.34. 76.2 and 121.92 m/s. The parameters that were 
used to characterize the wakes generated were again the velocity deficit and 
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the wake width. As noted in Section 5 the flow across the screens is not eni- 
form; in general, there exists either a diffusion-or contraction of the flow 
due to the fact that it is unrestrained. To accountfor this effect the ve- 
locity deficit and wake width ratios are plotted as a fuiction of contraction 
ratio and the value of the ratio taken at a contraction ratio of unity. An 
example of this process is shown in Figure 83 for the .95an thick honeyconb. 

The characteristic ratios of wakes for the various screens are shown in Figures 
84 and 85 with the screens listed along the abscissa in increasing K. In 
general, the higher the screen resistance, the higher the wake deficit sup- 
pression as predicted by the transfer fuictions shown in (29). In addition, 
the wake Inversion predicted for high resistance screens was observed in the 
41* open area perforated pUte configuration. The wake width was reasonably 
constant, and equal to unity as predicted by the model, although large varia- 
tions in this quantity were observed in the 41* open area perforated plate 
configuration. The accurate determination of the downstreai. wake width was 
parti cu-1 3 t 1 y difficult in this case a 

The model presented in the interim report (29) due to Taylor and Batchelor (38) 
is 
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Wake ” = 1 63 

Width 

Ratio 


and was used to produce the predicted values of the ratios as shown in Figures 
84 and 86. The measured value of the pressure drop coefficient, K, Figure 
and theoretical values of the flow angle ratio. « , (See Equations 66 and 67 
below) were used in these expressions. The flow angle ratio is the ratio 
between the outflow angle from a screen to the Inflow angle, bot e ng re 
enced to the normal to the plane of the screen. 
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Examining first the comparison of measured and predicted velocity deficit ra- 
tios, It Is apparent that the trend is well reproduced by the model. In par- 
ticular the contention that Taylor and Batchelor's gauze theory Is applicable 
to honeycomb and perforated plate made In (29) receives support from these 
data. In terms of the absolute level of the velocity deficit ratio, the pre- 
dictions are 1n reasonable agreement with the data. 



Figure 83 Variation of Veicx;itg Deficit Ratio Across 0.95 cm ffoneycoiAb With 
Contraction Ratio 
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is thus substantiated in this aspect ratio range. 
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verse velocity conponent over the vortex. This 7 »f ‘1« ‘rans- 
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change In the transverse velocity ^ 

the perforated plate, Indicating as mentioned l! ! !, 
diffusion process Is accelerated by the presence f fhT ’ 
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py the theory presented In the Interim report 
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The transfer function described there was 
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According to this model the suppression of transverse velocity component 
s solely due to the turning characteristic of the screen, which in turn is 
dependent on the pressure drop coefficient of the perforated plate 


•Jt = 1.1 (1 + 
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account for the diffusion process (vortex broadening) across the perforate, 
plates in order that the transfer function may be assessed, the upstream and 
downstream aaimuthal velocity fields are plotted as a function of transverse 

8 0 the 66* open area perforated plate. Subsequently the data is fUted bj 

minim zing the squares of the differences, to a potential vortex pro.'lle The 
profile used was ® 
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From this fit, an estimation of the circulation around the vortices upstrea, 

and downstream of the perforated plate is obtained. Now if the velocity field 
Is corrected for any bias— /3, then 
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Hence the ratio of the normalized azimuthal velocity fields across the per- 

slTand o‘'J?!T T “ '' «•- 

sored and predicted values of this velocity ratio is shown in Figure 87 This 

comparison is shown plotted against the mean velocity ratio across the 'rer- 
forated plates. ^ 
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Figure 86 The Vortex Trensverse Velocity Field Bach Side of the 66% Open Area 
Perforated Plate 
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Figure 87 Comparison of Measured and Predicted Transvetsh Velocity Field 
Ratios in a Vortex Across Perforated Plate 
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C<„sid.rin9 the difficulty in obteinina the true-value of the core radios f^ 
:e d7a e -pari son is enco^raain, and provides sore support for the « e,. 
mldel is thus sufficiently accurate for use in the ICS desi^ systen t 
ranains to provide sono estimation of the core radius change across per- 

f orated plate. 

The change in the core radius that is observed is shown graphically F'^ur® 
88 The diffusion of the vorter, not unexpectedly, is higher, the grea er e 
re;ista„ce of the perforated plate, and s.,.e indicati^ is 
curve fit of the degree to which the vortex core would be modified by 

trary porforatod plate. 

The transfer function of Equation 65 does not describe all 

Ten a vortex is convected through perforated plate. It • ' 

jmction with knowledge of the change in core radius, an estimate 

mutual velocity field change. 
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Associated with the transverse velocity field of the vortex Is an axial velo- 
city distortion. It Is cusp-like In shape and, as noted In Section 5, when the 
vortex Is convected through honeycomb, the deficit Is Increased (see Figure 
30). The cusp-like nature of this deficit Is most apparent after the honeycomb. 

The axial flow In the core of a vortex has been studied by, among others, 
Newnan (39) and Batchelor (40). Various experiments, see Ramsey (41), have 
produced conflicting descriptions of the direction of this flow. The various 
results described by Ramsey can be explained qualitatively by considering 
Satchel or^s work (40). A form of Batchelor's expression for the axial velocity 
field In a vortex is shown in Equation 69. 
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The axial velocity profile In the vortex Is a function of the static press^e 
distribution , the azimuthal component of the dynamic head 

and the total pressure loss distribution AH due to passage of air over the 
generating body. Depending on the relative magnitude of these contributions, 
the axial profile In the vortex may be a deficit (wake), a surplus (jet) or 
even a combination of the two. In the present experiments the axial velocity 
profile was a deficit. Indicating that the total pressure loss from the flow 
in the vicinity of the airfoil was the dominant term In Equation 69 since 
Batchelor has shown that the sun of the first two terms In Equation 69 is 
positive. 

Equation 69 Is now used as a starting point to determine the effect of a 
honeycomb on the vortex. The action of the honeycomb on this vortex Is to to- 
tally suppress the circumferential velocity field and thus render the radial 
static pressure gradient unsupportable. The first two terms In Equation 69 are 
thus ronoved and the axial velocity deficit Is Increased accordingly. 
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The expression for the ratio of the maximum velocity deficits across the 
honeycomb is developed in Appendix C and is 

max _ 1 . ^2A max 

AU, - X + 

“ lA max lA max 70 

Predictions using this expression are shown in Figure 89 together with the 
measured values of this ratio. The scatter of the measurements is possibly due 
to the error involved in estimating velocities at the critical points in the 
vortex. The shape of the axial velocity profile as predicted by the analysis 
in Appendix C is changed by the honeyccxnb (if the pressure loss across the 
honeycomb is uniform), but it is difficult to ascertain if the measured pro- 
files support this (see for example Figure 30), The axial velocity field is 
also modified on convection through perforated olate. 


^^lAmax 



figure 89 Vortex Axial Velocity Defi 'it Ratio Across Honeycomb 
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As recoNsd previously 1e Section 5. the 66* open erea perforated plate caused 
sn increase In the axial velocity deficit while the higher solidity (41* open 
area) plate suppressed the deficit. The measured values of the maxlmun axial 
velocity deficit ratio across perforated plate are shown In Figure gO. The 
theoretical predictions shown In this figure are derived using the model de- 
veloped In Appendix C for perforated plate, namely 
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l.A max 
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The data from the 66* open area perforated plate experiment Is reasonably well 
predicted by this model, however, the ratio from the 41* open area cate je of 
significantly lower value than the prediction. This difference Is probably due. 
to the higher resistance of this screen, which may cause the assiiaptlons of 
the.model to be v1olated« 

The perforated plate has a dispersive effect on the vortex as 
action on the azimuthal velocity field and the axial velocity f eld “ 

f«l the effect of this dispersion. The presence of the disperslm effect sug- 
gests that the pressure loss across the perforated plate Is not (as " 

L derivations In Appendix C) constant at every point In the vortex. It 
apparently vary radially such that (at least for the high 
plate) variations In the pressure loss dominate the velocity deficit ratio 
(Equation A.10). This effect seems to suppress the axial flow distor on me 
than predicted. Figure 90 and thus t-le theory provides an estimate of 
imun value of the downstream axial distortion. 

Frmn the above experimental results. It Is clear that the transverse velocity 
field in a vortex can be readily destroyed by use of honeycmnb. Tddt. 
vortices are likely to exist In the vicinity of any test stand 
include honeyconb. In addition, because the axial velocity distortion In t e 
vortex can be Increased by passage through the honeycomb. It Is Important 
put the honeyconb as far away from the fan Inlet as possible. This then per- 
mits slcy-lflcant flow contraction downstremn of the honeycomb to the fan face 
and as shown In the previous section 6.2.1.. the contraction will eliminate 
the surviving axial velocity distortion. 
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Figure 90 Vortex Axial Velocitg Deficit Ratio Across Perforated Plate 
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sumary, the models for the screens to be used in the ICS design systan are: 


1) Axi al Velocity Distortions 
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2) Vortex Vel oci ty Distortion 
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6.3 Models for the Acoustic Transmission Loss and Directivity Effects of 
Scr eens 


In the Interim Phase II Report (29) models were presented that enabled the 
effect of the ICS on sound propagation to be determined. Specifically, expres- 
sions for the transmission loss of screen types (honeycomb and perforated 
plate) were shown... Also, a model of the effect of a corner made of honeycoab 
panels was developed and shown to agree closely with measured data. In this 
section the transmission loss model is assessed using the data gathered during 
experiments described in Section 5.2. The corner diffraction model was applied 
to a honeycomb corner in the Interim Phase II Report (29), and is here applied 
to structural manbers and the results are applied to measurements of the 
directivity effect of solid structure* described in Section 5.2. Finally, the 
impact of these effects on the measuranents of farfield sound pressure due to 
an engine operating with an. ICS is discussed. 

6.3.1 The Acoustic Transmission Loss of Honeycomb Panels and Perforated Plate 

The transmission loss of various screen materials was determined as described 
in Section 5.2. As noted there the large standard deviation of the measure- 
ments coupled with the small transmission loss of the material indicated that 
there is no significant transmission loss for many of the materials and flow 
conditions examined. 


The model for the transmission loss of a sound wave normally incident to the 
honeycomb. 
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predicts an effectively zero attenuation at all frequencies. The measured 
transmission loss mean for the .95 cm thick honeycomb at 12.19 m/s is shown in 
Figure 91. The standard deviation bounds are also shown and bracket the pre- 
dicted value of zero. Tl)e large standard deviation coupled with the measured 
mean being close to zero confinn the predictions of the model , i.e., that the 
transmission loss is small. The results from this configuration are typical of 
those from the 1.27 on honeyconb experiment. In Section 5.2.1, it was shown 
that there is an increase in mean transmission loss with flow speed for the 
thickest (2.54 cm) honeyccmb. This is indicative of a non- zero resistive com- 
ponent in the honeycomb impedance. This characteristic is not accounted for by 
the honeyconb transmission loss model. Now since it is expected that the 
acoustic resistance is a function of the pressure drop coefficient, K, of the 
honeycomb, the range of applicability of the transmission loss model can be 
defined by K. From the measurements of K (Section 5), the thinnest honeycombs 
have pressure drop coefficients of .41 and .54. Thus a honeycomb having a 
pressure drop coefficient less than about .5 operating in a flow of less than 
12.19 m/s will attenuate a sound field according to the model. Equation 79. 
This model could underpredict transmiss'on losses at higher speeds awLpres- 
sure drop coefficients. Equation 79 can be manipulated to provide a guide for 
honeycomb design that minimizes the transmission loss. If it is desired that 
the attenuation due to honeycomb be less than 1 dB, It is necessary that 

sin n sin k ^ V 

< 1.03 80 

cos (' 


For - 1 , the worst case, this inequality shows that the incidence 

angle for the sound field can be as high as about 50° without resulting in 
an attenuation in excess of 1 dB. The implication for the ICS design system is 
that as far as honeycomb transmission loss is concerned, a Sf^erical design 
will satisfy the condition of Equation 80 provided that the ICS radius is 
greater than about twice the fan dianeter. 
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The model Indicates that It is possible for some perforated plates to signifi 
cantly attenuate sound at the higher frequencies. This can be prevented by 
ensuring that the material satisfies the criterion. 
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which establishes a maxirawn trahsmisslon loss of 1 dB In a frequency band 
bounded by fj and tj. Any deviation of the Incident wave from nomal will 
nrnduce a lower transmission loss In the perforated plate. 

The difference In behavior of the two perforated plates In frequency Is not 
explainable In terms of the model. The higher transmission loss of the lower 
solidity plate is also confusing. However it should be noted that neither the 
screen of (15), (52* open area) nor the perforated plates of (22) exhibited 
any observable transmission loss. 


While the honeycomb itself is acoustically transparent a corner formed by ad- 
Jacent panels or structure can affect the directivity pattern in the far field 


6.3.2 Directivity Changes Due to a Honeycomb Corner and Structural 
Members 


Measureirents made of the sound field behind a honeycomb corner were shown In 
the interim Phase II Report (29). A theoretical model of the phenmenon was 
also presented there and good agreanent was observed between measurement and 
predictions. The experiment described In Section 5 also Included the 
definition of the directivity changes due to structural members. The struc- 
tures examined were a square rod and a circular rod. The diffraction mode 
used to predict the honeycomb corner effect on the sound field was adapted to 
predict the effect of the rods on the field. Comparisons between measured val- 
ues of the Change in sound pressure level due to the presence of the rods are 
shown together with predictions in Figures 94 through 99. Both pattern shape 
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and absolute level of the directivity change are well predicted. The predic- 
tion Is at most about 1 dB different from the measured data. It Is Interesting 
that the shape of the rod cross-section does not affect the gross characteris- 
tics of the far field. The Important parameter Is the dimension of the ob- 
struction normal to the Incident wave number vector. The model then developed 
In the Interim Phase II Report (29) gives good estimates of directivity chan- 
ges due to both corners In honeycomb panels and structural members. 

There are two Important differences between the Idealized model and the actual 
situation at the test stand. Firstly, the source Is not a point (as assuned In 
the model) and secondly the source does not radiate at a single frequency but 
over a wide band of frequencies. The actual source noise spectrua can, how- 
ever, be considered to be made up of a combination of discrete tones and broad- 
band noise. 

The radiation field of the discrete tone component can be considered to be due 
to an array of point sources distributed over the engine Inlet. Using the 
model and a distribution of sources along a line, the deviation from the free 
field radiation pattern due to the presence of a corner was computed. While 
the line source cannot be considered a simulation of the actual distributed 
source of the engine, the resulting directivity patterns of the two source 
types are qualitatively similar. In consequence, the distributed line source 
(and the two source models examined below) Is useful in defining trends. The 
source amplitudes of the line source were uniform and an arbitrary source 
phase distribution was considered. The honeycomb corner simulated one found 
on the P&WA developed ICS. The schematic of this simulation Is shown In Figure 
100. The resulting directl vity modification (I.e., the difference between the 
directivity patterns with and without the honeycomb corner) 1s shown In 
Figures 101 and 102 for frequencies of 1850 Hz and 3700 Hz, respectively. 
(These frequencies correspond to JT9D blade passing tone and harmonic fre- 
quencies at an approach condition.) For the purpose of canparison the direc- 
ts vity modifications due to the same corner with a single point source are 
also shown In these figures. The comparison Indicates that the radiation field 
of the more realistic distributed source Is disturbed more by the corner than 
that of a point source. 
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figure 94 


The Directivity Modification Due to a Square Bar 
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rl«ure 95 Tte Dlrectivita Modification Due to a Square Bar 
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Figure 98 The Directivity Modification Due to a Round Bar 
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Figure 99 The Directivity Modification Due to a Round Bar 
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Figure 102 The Directivity Modification of a Honeycomb Corner at 3700 Hz 


In contrast, if the source phases are allowed to vary, a different picture 
anerges. Consider, for illustrative purposes, a distributed source comprised 
of two sources. Figure 103. If the phase of one of these sources is varied 
with respect to the other and averages performed over the radiation patterns, 
an average directivity modification pattern can be computed (This changing of 
source phase relationships occurs In practice with variation of fan speed, and 
the averaging process over phase, described above, has a counterpart in aver- 
aging over engine speed in the test situation). This computation has been per- 
formed for several cases. 
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Figure X03 Diagram of Discrete Frequency Distribution Source Simulation 


In Figure 104 the average directivity change has been computed for source var- 
iations over 1°, 2° and 10*^ at a frequency of 1850 Hz due to a corner at 
67.5®. The averaging process produces a noticeable "washing out" of the cor- 
ner effects for very small source phase variations. A similar plot is shown in 
Figure 105 for 3700 Hz and the same observation can be made at this frequency. 
The power of the "washing out" of corner effects by small source phase changes 
is apparent from this simple simulation. For the more accurate approximation 
of many sources distributed over the inlet disc,, it is expected that small 
variations in the phase of the individual sources will suppress any corner 
effects even further. The impact of this observation on the ICS design system 
will be assessed in the following subsection 6.3.3. 
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Figure 104 The Effect of Averaging Over Source Phase Differences on Corner 
Directivity Correction (2 Sources./ 
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The ecouetlc eiellarity pereeeUrs^^. \ and i are e^wh f'^uras 
Z, 102. 104. and 106. evaluated for the 0T9D approach condi t.on^and a sta ic 
test configuration. For the OT160 these values of the paraneters^ (no a 
,«d microphone radius). -H (nomalited ICS corner radius) 
honeycomb thickness) are t?^ical. In particular the honeycomb t ckness 
quite realistic. The only variation from the JW case then is in 
values Of this paraaeter for the 0T15D *t - J»ach -^^i- -d at 
Dassing frequency and its hamomc are 7.14 and 14.28 respect y 9 

1 *t y the directivity modifications to the radiation field, in general. 
:lcfeale in magnitude, the higher the value of ^all other P---- e,n^ 
held constant. Consequently since the JT150 values of 

corresponding values for the J190. all of the previous ccmmen s addressed 
the JT90 configuration are applicable to the 0T150 configuration. 

The second conponent of the fan source spectrum is the broadband cosponent. 
for a rando. so^ce. the amplitude and phase of the source at > 9'«" ^ 

quency are randan variables. The consequence of this is that, a 9 
frec^ency. the directivity change due to the corner «ill tend to be .washed out 
byThe rLonness of the soiree in a manner similar to that discussed pre- 
viously for discrete tones. Another factor in determining broadband correc- 
tions is that they are defined on a one-third octave band basis. The thir 
octave band correction is conceptually determined by integrating the intensity 
fields over each frequency range with and without » cor-r in pl.c^ and ev - 
uat -0 the ratio of the resulting one- third octave band intensities at each 
Z J ?lTe stmple simulation of Figure 103 was used to observe rffect on 
th! directivity change due to simation over frei,iency. In Figure IM a com- 

‘p'rirrs lad: betwfen tim directivity modification of the P-- - 

^ j- V- nhaco at a sinale frequency of 4000 and tnac 

sence of source radiating in phase at a singie j 

Obtained when the soirees are radiating in phase over the 

tave band (3665 Hz-4488 Hz). (The one-third octave 

tamed by approximating the integral over frequency J 

pari son shows the reat oaount of suppression of corner effects that occurs 

when the modification of a broadband radiation field is examined on a one- 

third octave band basis. 


131 


The broadband directi vity modification due to the honeycomb corner is minimal 
because of source randomness in phase and amplitude and the integration over 
one- thi rd - oct ave bands . 

The consequences of these conclusions for the ICS design system and static 
testing procedures are discussed in the following section (6.3.3), 

6.3.3 The Impact of the Acoustical Transmission Effects of an ICS on 
ICS Design and Static Testing Procedures 

The transmission loss of typical possible ICS component materials (honeycomb 
and perforated plate) is negligible provided the design criterion presented in 
Section 6.3.1 are met. These criteria are for panels having low through-flow 
speeds. Corners between adjacent panels and ICS support structure can poten- 
tially modify the far radiation field. 

It has been observed that an engine fitted with an ICS operates at a given fan 
speed with little, if any, variation from that speed. The stability of the 
engine operation is greatly enhanced by an ICS. Associated with this stability 
of engine operation is an increase in the stability of the discrete tone 
source. This is manifested in the far field by a much more peaked discrete 
tone field. However, even in this more stable operating condition, the tone 
level at a given angle varies greatly because the increased tone level grad- 
ients in the far field make the tone level very susceptable to any small var- 
iations in the source and/or the radiation field. The model is used to illus- 
trate this point as shown in Figure 107. Here a part of the free radiation 
field of two sources radiating in phase is compared with that obtained when a 
0.3° phase difference exists between the sources. At some angles the level 
changes radically (lOdB) while the general shape is unchanged. Inlet kulite 
data gathered during the Joint Noise Reductioti Program (JNRP) also exhibit 
some unsteadiness in tone level with time, with an ICS in place, although the 
level of unsteadiness is lower than with no ICS in place. It is possible that 
this unsteadiness encountered when testing an engine with an ICS is sufficient 
to eliminate any corner effects on the radiation field, however, a more con- 
trolled and certain suppression is desirable. 
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Figure 106 


Figure 107 
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Variation of the Free Field Radiation Pattern of 2 Sources Due to 
a Small Phase Change 



AS demonstrsted In Section 6.3.2 n.od1ficotions to the discrete tone redUtlon 
field by a corner can be suppressed by averaging over source phase variations. 
The phase (and amplitude) relationships In the distributed source at the en- 
gine inlet may be varied by changing the fan speed. This has 
for small mordents In fan speed over a small speed range, on a Tgo t ed 
with an almost cornerless ICS (developed by The Boeing Company). Correlation 
techniques were used to detemlne the change In fan speed necessary to des roy 
a^rsiliUnty between blade passing tone far field directivity patterns. 1 
was found that a fan speed variation of the order of .5* produced uncorrelated 

far field patterns. 

A similar analysis was performed on the patterns produced by the two source 
simulation examined In Section 6.3.2. A aero correlation between source ra - 
tion patterns was reached after a variation of the source phase difference 
2«. AS can be seen In Figures 104 and 105 the directivity modification to 
the radiation field by the corner when averaging over 2 is sigmflca.tly 
reduced to approximately a maximum of 1 dB. It Is therefore suggested that 
directivity modifications to the discrete tone radiation field due to e pre 
sence of the corners and structural members of an ICS can be eliminated by 
varying the fan speed continuously during data acquisition by an amount 
order of .5% of the mean fan speed. 

Directivity modifications to the broadband radiation field are probably neg- 
ligible due to 1) source phase and anplltude randomness (even during stable 
engine operation), and 2) the practice of quantifying broadband noise leve s 
ever one-third octave bands. It Is considered that either of these causes is 
sufficient to eliminate broadband directivity modifications due to pane 
corners and structural members. 

U Should be noted that the honeycomb corner used In the Illustrative simula- 
tions of section 6.2.3 has the characteristics of the PW developed ICS. The 
constructi on techni quo of thi s I CS resul ted 1 n d1 sconti nul t1 es caused by 
corner structure and .loints In honeycomb panels at these corners. Subsequent 


ICS Designs (e.g. those of The Boeing Company and NASA Lewis) minimized the 
severity of these corners by using Improved construction techniques. Conse- 
quently, the magnitudes of directivity modification due to a corner presented 
here, based on the P&WA ICS design, provide an Indication of the upper bound 
of such change. Moreover, In spite of the severity of this ICS corner, no ef- 
fect on the radiation field that could be attributed to the presence of a 
corner has been observed In P&WA experimental data, A possible explanation for 
this is the coarse measurement grid (every 10®) used to define the radiation 
field. From the patterns generated using the directi vity model , extreme direc- 
ti vlty modifications occur over narrow far field angle ranges and the proba- 
bility that a microphone In a coarse grid will appear at such an angle Is not 
great. Coupled with this, any unsteadiness In the ambient conditions would 
contribute to the reduction of corner effects. 

Therefore, It seems that even if severe discontinuities exist in the ICS con- 
struction, any directivity modifications resulting from them can be suppressed. 

However, as discussed in the Interim Phase II Report (29), the problem Is 
preferably eliminated at the design stage by minimizing all such discontinui- 
ties from the outset. The Boeing developed ICS Is an excellent example for 
full scale engines as Is the NASA Lewis produced ICS for the JT15D (23). 
Acoustically these designs are highly satisfactory. 

6.4 The ICS Design System 

In the Interim Phase II Report (29) the main elements of an ICS design were 
Identified. They are: 

0 The criterion for material selection. 

0 The ICS self generated distortion constraint. 

0 The ICS acoustic transmission loss and directivity constraints. 
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. nf this section Is to aenerete a design procedure, that Incorpor- 
rreirtl^ u™s . par.eters readily aeanahie to the 
signer. The models detailed In the previous sections fora 
ICS design system. 

a.4.1 The ICS Location and the Characteristic Dimension of the Screen Material 
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tour placed in the engine dmnlnated fl<« reg.on th 
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I'^/o oTs. tI would violate the constraints discussed at the beginning of 
this section. 
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consequently. It is suggested that the ICS be placed in the third regim and 
be designed so that the engine induced velocity is 1.5 ti^es the «axi« a^ 
bient wind speed in the test nindow. This ensures that air fl»s into the IC 
at all points and also that the pressure drop across the ICS is snail. Corre- 
spondingly any total pressure or strea»w1se velocity distortion do«nstre» ^ 
the ICS would be snail and readily suppressed by the flow contraction. I 
unlikely that the velocity distortion would generatf significant turbu ence. 
in view of the low distortions generated at such a surface in this region, 
actual shape is probably not of great importance as long as it is roughly 
hemispherical. This insensitivity to shape is borne out by BMI measurements 
made with the PWA. Boeing and NASA Lewis ICS-s installed, where there was no 
evidence of distortions introduced due to the shape of these ICS s. The radius 
of the ICS can thus be determined from the engine inlet velocity at its lowest 
operating condition and the maximim. anbient wind speed in the test window 
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The characteristic dimension of the screen material can new be detemined. 

The minim™ turbulence transverse scale to which the engine is sensitive will 
have a corresponding value at the ICS surface 

1/2 t 


( ^min ^ ICS 1 
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where the contraction ratio is that at the lowest engine speed and maximim 
rbient Wind speed. How in order for the ICS to act «"«ti«ly. Us charac- 
teristic dimension should be roughly an order of magnitude less than the m 
mim value of this scale at the ICS. This will ensure that turbulence scales 
likely to produce BPF tone noise will be suppressed by the screen. 
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Thus, 
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This dimension refers to the solid material of the screen as well as the hole 
diameter of both honeycomb and perforated plate or gauze. 


A flow chart of these calculation processes Is shown In Figure 108, 



Figure 108 Flo» Chart for Determining ICS Location and Detail Dimensions 
6»4.2 The Criterion for ICS Material Selection 

Two possible criteria for screen material selection were suggested In the In- 
terim Phase II Report (29). One required that the material be chosen such that. 
In sane sense, the Inflight fan face flow distortion field was simulated. The 
other criterion suggested there was that the screen material be chosen such 
that all fan face Inflow distortions generated upstream of, and by the ICS, be 
totally suppressed. Both of these criteria were confined to those elements of 
the distortion field outside of the nacelle boundary layer that would generate 
tone noise on interaction with the fan. 
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These criteria are, each for Its own reason, difficult to attain. The simula- 
tion criterion. In addition. If attainable would strictly require a different 
material for each operating condition while the total suppression of fan face 
turbulence Is probably Impossible. 

It Is therefore suggested that a more attainable criterion be that the ICS 
material be chosen such that, In some sense, the fan face Inflow distortion 
encountered statically be less than or equal to that expected In flight. The 
distortion field at the fan face has previously been considered In two parts a 
turbulent component and a steady (or quasi-steady) distortion. Following this 
division, the quantity that is best suited as a measure of the turbulence 
field at the fan face will now be examined in order to formulate the turbu- 
lence distortion criterion for ICS design; subsequently the steady (or quasi- 
steady) distortion component will be examined. 

The detailed characteristics of a turbulent field are provided by the spectrim 
tensor, a quantity most easily explained as the Fourier transform of the cor- 
relation tensor. The diagonal elements of this tensor may be interpreted as 
the distribution of the energy associated with each velocity conponent in wave 
nimber space. If it is assumed that the rotor- tucJxil ence source is primarily 
due to the fluctuating lift on the fan blades then the element of the spectrum 
tensor on which the source strength is most dependent is the upfyash element, 

>’nn 

This may be considered as a surface, infinite in extent. Now Pickett (42) and 
more recently Ganz (36) have implied that not all of this surface contributes 
significantly to the level of BPF tones (i.e. only the energy in certain wave 
number ranges contributes to BPF tones). This tone generating part of the 
spec- trim surface may be defined by the inequalities. 


^2L " ^^2 " ^2U 
'‘31, '‘J ■ '‘il- 
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It should be noted that the three dimensional spectrun is discussed here in 
terms of Cartesian wave nimbcrs, whereas in view of the inlet gecmetry, t*he 
polar versions are more appropriate; however, it is assjned that locally the 
field may be represented in the Cartesian fr<nie. There are then only a block 
of wavenimbers at which turbulent energy produces BPF tone noise. Figure 109. 
(In this figure only two wave nuaber conponents are represented, however the 
surfaces is in all three components). The characteristics of this block are 
dependent on engine size, fan blade nimber and engine speed. In terms of eddy 
sizes, this wavenumber block defines eddy geometries that can generate RPF 
tone noise upon interaction with the fan. A possible design criterion is, 
then, that the turbulence energy density in this wavenumber range be less 
statically than in flight. A simpler but approximate form of this criterion is 
that all three one dimensional PSDs based on the wavenimber range of interest 

be less statically than in flight in the appropriate wavenunber ranges. This 
would produce three designs in general. The differences in these ICS designs 

reflect the difference in shape between the flight and static fan face upwash 
spectra. If the shapes are similar the three screen designs would be similar 
and the design system could be based on one PSD alone. 





Interact ion Kith a Rotor 
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Assuming this to be the case, the one dimensional PSD ‘I'l 

wave nunber range of interest will be used in the ICS design criterion. Thus, 

the criterion for turbulence suppress! m is 


>’n<>'2>S ^ fn"'2'p 


kj < k^* 


^2L ^2 ^2U 


k3L < k3 < k3^, 
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To determine the target value , the flight fan face upwash PSD, it 

will be assumed that the contraction in flight is negligible. Consequently, 
the fan face turbulence field will be assimed locally isotropic in the desired 
wavenumber range in the manner of the local flight atmospheric turbulence 
field. That field is best modelled on a one- dimensional PSD basis by the Von 
Karman spectrum (36). However, to evaluate ^^^^2^? 

wavenumber ranges of inequality in Equation 86 requires knowledge of the three 
dimensional spectrum ^ this manipulation the spectrum tensor that 

produces the isotropic Dryden spectrum (36) is more suitable. The upwash ele- 
ment of this tensor is 

2 2 

'nn'!S' = '’ll >o + ''22 COS Y^+ CROSS TERM 87 


0*" r (k) = 
nn 


N. 




j^k^ cos Yo + k2^ sin^ ''"o ^ j ^ 
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This is derived in Appendix D. 
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Tho constants and are chosen such that In the wavenunber range of 
interest the Dryden spectrim approximates the atmospheric Von Kaman spectrun. 
The determination of N^, and Vp is shown in (36). The relationship between N, y 
and the variance and integral length scale of the turbulence are discussed in 
Appendi x 0. 


From desired one dimensional PSOFj^ik^lp can be detenained from, 


k k * 
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With the aid of suitable approximations this integration can be performed and 
is shown in Appendix 0. The result is 

(’^30 ■ '^ 31 , • ' - - - 


|'4k2^ sin*" 

h k 

3U ^3L. 

2 



^ ^3U '^3L , 

2 

3 


The integration limits and k^* are related to the range of 

turbulence scales over which, upon interaction with the fan, blade passing 
tones are generated. Thus, if the turbulence scale ranges to which an engine 
is sensitive and the atmospheric turbulence parjmeters in flight, Np and Yp , 
are known, the target flight fan face upwash PSD i^ri(k 2 )p can be determined. 

It is now necessary to formulate the second part of the inequality 86, the fan 
face upwash PSD ^^ 2 ^‘^cncountered statically with an ICS. To do this 
it is helpful to have a view of *^he operations and intermediate states that 
the atmospheric turbulence go*" through, statically, on its way to the fan 
face. These are shown schematically in Figure 110. As in flight, the atmos- 
pheric spectrum tensor is approximated by the isotropic tensor ^ in the 
wavenunber range of interest. This tensor is defined upon the appropriate 
choice of N^, and 'i s (See Appendix 0). Subsequently this turbulence undergoes 
contraction to the ICS. 
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Q DENOTES A STAVE 
□ MENOTES AN OPERATION 


Figure 110 


Flm Chart of the Processes to Which 
Subjected in Static Engine Operation 


Atmospheric Turbulence 


is 


The nature of contraction fr<™ the atmosphere repul res seme discussion In view 

of the possibility of an Infinite contraction ratio. Firstly, It is necessary 

to distinguish between turbulence present In the atmosphere independent of 

engine operation and that which relies upon the engine for Its existence It 

Should be recalled that this section Is concerned solely with the former'type 

On running an engine, then, this turbulence Is Ingested into the fan and for a 

given engine condition and ICS configuration the contraction ratio from the 

atmosphere to the ICS Is deteimined by the local mean- atmospheric wind speed 

If an engine is run in the FAR-36 test window, this contraction ratio can vary 

fron ■■ to small values -(around 3). It Is consequently necessary to choose 

the most appropriate value of the Initial contraction ratio for the ICS design 
system. ^ 


For the Identification of this contraction ratio consider the approximate one 
dimensional post-contraction PSDs presented in (36). An approximation Is per- 
formed there by Ganz fort • i , where. ^ 

This approximation Is satisfied in the static test window. The post contrac- 
tion PSD In terms of the pre-contraction wavenimber Is given there as 


14 .? 






where » and p are functions of y ~ . 

Note that the shape Is determined by and 1’ and thus by alone. 

Now the pre- cor.tr action wavenuiiber, , where k^j^is the post 

^ L _ '^1 . 

contraction wavenunber and ’ I'ua . which from Ganz's model of the Ingested 
atmosphere (36) Is a function of the contraction ratio by virtue of the samp- 
ling of the atmosphere by the capture streamtube, so that 
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Thus 




k 


2n 


Uj!') U 
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is independent of contraction ratio. Correspondingly the shapes of 
the post contraction PSOs in the transverse wavemntbers are independent of 
contraction ratio. In addition, a small contraction is sufficient to transfer 
most of the energy in the streaniwise wavenumber PSDs into the range where the 
blade passing tone level is independent of axial length scale. In short the 
tone generation capacity of the incoming turbulence is, above a low contrac- 
tion ratio, independent of the distortion of the turbulence field and is 
greater or lesser solely according to the post- contraction upwash energy. The 
initial contraction ratio must therefore be chosen to account for the highest 
level of upwash energy encountered at the ICS in the test window. 


The mean square value of the upwash velocity component at the ICS is 
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u , 




t u , 
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which may be written in tenns of the moan square value of the atmospheric tur 
bulence velocity (isotrofic initial conditions assuned) 
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where i'| and "2 are the mean square velocity component ratios across the 
Initial contraction. Now using Prandtl's approximations 



" '1 
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andGanz*s (36) expression for the initial variance of the atmospheric turbu- 
1 ance field sampled by the capture streamtube, 
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in Equation 94 yields. 
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Thus, the post- contract! on upwash kinetic energy decreases monotonically with 
contraction ratio. In consequence, the contraction ratio in the test window 
that produces the most severe turbulence field at the ICS is the lowest value 
in that window. Knowing this contraction ratio, together with the initial at- 
mospheric turbulence conditions, the model sumnarized in Section 6.1 can be 
used to determine the turbulence characteristics at the station upstream of 
the ICS, provided that the viscous dissipation factor can be quantified. 


Since the lowest initial contraction ratio in the test window occurs at the 
design point, i.e., = 1.5, the contraction distance will be small compared 

to the distance from the origin of the turbulence. Consequently viscous decay 
effects are negligible. 


At the entrance to the ICS then, from Equations 28 and 29 
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the initial contraction r 


probably be less than 6. 


since oo 

icales at the ICS from Equations . i 


correspondingly, the tntecyal lengths, 

33 , 34 and 35 are 
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The ICS material can be described prersure' drop coefficient), 

paraaeters. n (tl. fl»- angle rat,o and K^,^(th ^P 

These paraaeters have been discus characteristics, 

,^oe„ ™odel Of section h.l „™,ediately behind the ICS (negiect- 

the turbulent flow-fie i c ai square values of e 

,„,self ,neratedtur^lence)»ay he found. 

turbulent velocity coaponents there are 
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. , lenqth scales are unchanged and 
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,S then contracted to the fan face. 
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At this point It is necessary to make some comments concerning the inflow con- 
trol structure mater i al. If the inflow control structure is composed solely of 
honeycomb^ the transverse vel ocity components will be effect! vely suppressed 
while the streamwise turbulent velocity components will convect through the 
honeycomb largely unaffected. In consequence the turbulent energy distribution 
will be strongly weighted in favor of the streamwise velocity component, inv 
mediately downstream of the inflov control structure. However, due to a) the 
tendency to isotropy in turbulent flows and b) the amplification of the trans- 
verse velocity components of turbulent flows under contraction, there is prob- 
ably a rapid increase in the amplitude of the transverse turbulent velocity 
component downstream of the inflow control structure. Now if the extraneous 
noise generating inflow distortion is a steady or quasi steady vortex then it 
is possible that a honeycomb may be sufficient, however it is probably not 
possible to a) make this assessment of the inflow distortion field aprior or 
b) if it can be made, ensure that this assessment remains valid under all 
conditions of use of the ICS— In t his case, the use of a resistive elenent in 
the ICS would provide a conservative design. It is therefore, in general, 

necessary and/or conservative for an ICS to suppress both transverse and 

streamwise distortion vel ocity components i.e. the ICS should include both 
resistive and torque absorbing elenents. 

N(m the contraction model of Section 6 requires isotropic turbulence as an 
initial condition and so to use this model, the design condition of isotropy 
is imposed on the ICS. Thus, while many ICS designs could satisfy the design 
criterion of producing an upwash field at the fan face in static operation 
less than or equal to that encountered in flight, the one that produces equi- 
partition of the turbulent kinetic energy downstream of— the ICS is chosen. In 
general, this will require a resistive element (gauze or perforated plate) and 
a torque absorbing elanent (honeycomb). Equi partition of the turbulent kinetic 
energy requires that downstream of the ICS 



If it is assumed that the two elements act Independently and in series, then 
from Equation 101 
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In view of the small initial contraction the integral length scale relation- 
ships will exhibit roughly isotropic characteristics and a scale representa- 
tive of this field is 
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For the purposes of the final flow contraction the equi partitioned turbulence 
field downstream of the ICS Is then, approximated by an isotropic field of 
characteristic length scale =" f -22 As noted previously the upwash field 
resulting from the contraction of an isotropic field is dominated by the 
transverse (azimuthal) velocity component in the contraction range encountered 
on a test stand, hence 
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and *^2^^2^s can be obtained from Ganz's extension of the Ribner-Tucker 
theory (33), thus 
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The ICS subscript refers to the post-screen location and k, is the local post- 
final contraction wave niitiber. 


With the formulation of‘n 2' sit is now possible to expand the design cri- 
terion of Equation 86 with the help of Equations 90 and 105 so that the cri- 
teri on becomes 
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This criterion together with the honeyconib and perforated plate (or gauze) 
characteristics and the equipartttion constraint of Equation 102 constitute a 
0 our equations in the four unknowns “p ' f K_ , k 

(the characteristics of the teo screen elanents). The honeycnb characteristic 
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Figure 111 The Elements of the ICS Material Selection Criterion 




















The other element of the external field that must be considered is the steady 
or quasi steady velocity distortion. These can be in type either streamwise 
velocity deficits or vortices, Their origins can be surfaces over which the 
air travels on its way to the engine or mid-air stagnation points. Stand 
structure can be a source of both wakes and vortices, since variations in flow 
incidence angle over the structure can cause local separation over the body 
and subsequent vortex roll-up. This phenomenon has been observed in bluff 
bodies, for instance (43). The variation in loading over the support structure 
leads to the formation of potential vortex distributions. In addition, the 
ground plane can produce an ingested vortex and also the lesser known vortex 
originating at a mid-air stagnation point (44). Protuberances on the external 
nacelle have also been observed as sources of steady (quasi steady) inflow 
distortions. 

The quantification of the initial conditions of these flow disturbances is 
difficult in view of the conplexity of test stand geometries and more work 
needs to be performed in order to define them. 

As with the atmospheric turbulence, a criterion for suppression of this type 
of disturbance must first be established. The steady inflow distortion in 
flight is ill-defined and is due if it exists for different reasons (e.g. air- 
craft angle of attack) from those encountered statically. It is therefore sug- 
gested that the amplitude velocity components of the steady inflow distortion 
field be reduced below the RMS value of the atmospheric turbulence field en- 
countered in flight, or 
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where A indicates the anplltude of the distortion. The right hand side of 
these inequalities is readily determined as in the first part of this section. 
For the determination of the left hand side of the inequalities 110 and 111, 
it will be assumed that the streamwise velocity distortion is suppressed by 
the ICS Independently of the suppression of the transverse velocity distortion. 


If the tMtia) condition of the streamwtoe distortion is ,n . 
initial contraction, usino the mnrfoi f after the 

at the ICS entrance wiiTL ‘-'a deficit 
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and assuaiing the heneyconb and perforated plate (oaeze) an i„ • 

deficit downstream of the IfS is « <9aete) act m series, the 
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9 n as in the case for atmospheric turbulence inaestion th ^ 

wise velocity distortion will a.u ^ stream- 

y aisiortion will occur at the maximun anbient veloHtty ,-n 

test window, unless the wake generating structure lies ^.r ! 

by the engine inflow (which U unlikel^ 

The streamwise distortion criterion in Equation no is then 
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153 


On pursuing a similar course with the transverse (vortical) distortion, the 
inequality of Equation 111 be cones 
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Note here that the design criterion is independent of contraction ratio due to 
the invariance of a vortex under contraction. These two inequalities, together 
with the honeyconb and gauze (perforated plate) characteristics, fom a set of 
four relationships in four unknowns which can be solved. A schematic of this 
aspect of the design systen is shown in Figure 113 and the step by step pro- 
cedure is given in Section 6.4.4. 

6.4.3 Other Desi gn Aspects of an ICS 

various observations on the general features that should be incorporated into 
an ICS have been made in this report and (29). In sunmary these are: 

0 Discontinuities whether of structure or profile should be minimized in 

order to minimize sel f- generated distortions and directi vity modifications 
to the radi ati on f iel d. 

0 The perforated plate (qaiite) transmission loss criterion of Eqaation 82 
stmuld be applied to determine the allowable plate (gauze) thickness and 
hole diameter combinations. 

0 The construction of the ICS aft of the inlet plane is an open question. 
Both baffles and ICS screening material have been used in this region. The 
nacelle boundary layer is undoubtedly affected by the conditions here. 
However, evidence to date gathered under the JNRP has indicated that there 
is no significant difference in the noise fields of either configuration. 
In contrast, the work perfonned in (21) showed that the noise field was 
quite strongly dependent on the flw field in the boundary layer. The 



Figure 113 Flou Chart for Determining ICS Material (Steady Distortion 
Criterion J 
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sensitivity of the noise field to the nacelle boundary layer undoubtedly 
varies from fan to fan. It would be expected that an ICS open to the rear 
would allow more accurate noise measurements to be made at higher <"ar 
field angles. Nevertheless, the whole question of simulation of the flight 
nacelle boundary layer, in static operation, needs more Investigation. 

0 The use of two ICS materials In conjunction leads to the possibility of 
flow blockage when holes do not match up. This problem can be avoided by 
separating the two materials or alternately ensuring that the mesh size of 
one of the materials Is significantly greater than the other. This latter 
solution is quite acceptable in view of the degree of latitude allowed to 
the detailed dimensions of the ICS material (Equation 85). 

0 The relative location of ICS materials, I.e., whether perforated plate 
should be placed upstream or downstream of the honeycomb Is probably un- 
important as regards the suppression of the noise generating elements of 
the inflow field. It may be that because of construction requirements one 
material may produce greater downstream distortion than the other. In this 
case, an indication of the best configuration is provided. 


6.4,4 The Step By Step ICS Design Procedure 

1. Determine inlet velocity at lowest operating condition - 

2. Determine maximum ambient wind velocity in test window - 

3. Conpute design speed at ICS - 

4. Determine Inlet radius - 

5. Compute nominal radius of ICS - *^ics 
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6 . 


7. 


Determine blade number - 

Compute minimun sensitive transverse scale 
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8. Compute final contraction ratio (ICS to inlet) - IP 

II 




IV ■ II 


o 




9. Canpute minimum sensitive transverse scale at ICS - ( \n,in } tcs 


( ^min ) TCI 


IF 
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10. Compute characteristic detail dimension of ICS material - 

IC.S 
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ICS 
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\nin ) ICS 
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11, Compute Reynolds Number of honeycomb cell - Re 


Re 


u 


ICS ICS 


Atmospheric Turbulence Design 

12. Determine engine height - ?: , roughness scale . static reference 
height (engine height) z , static reference wind speed (max, wind 
speed in test window) 

Rr*F 

13, Compute RMS value of turbulent velocity in capture stream tube during 
static operation - 




14. Compute static contraction ratio - ' 

I s 

u 

o 

^ 1 s " if 
A 

(Note initial contraction ratio, from atmosphere to ICS, ' , , = 1.5 
by des i gn ) * ^ 
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15 . 


Ccmpute characteristic length scale of turbulence in capture 
streamtube during static operation - 

ii i’ ^ 

’ % 'is 


16. Determine, for flight operation, engine height £ , roughness 

scale , flight reference height - . flight reference wind 

speed - (I 

ur:p 


17. Compute RMS value of turbulent veXocity in captuv'e stream tube in 
flight - / ■ 

■ ^ 
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18. Compute characteristic length scale of turbulence in capture 


streamtube in flight - 
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Compute limits of sensitive wavenumber ranges - k,*/ k, ^ , k., 
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Compute in flight turbulence field characteristics - Nj, 
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Determine blade relative inflow angle jt tip - 

Compute characteristic length scale of turbulence behind ICS 
4 ' 'AS 
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23. Compute the parameters p and s 
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24. Compute maximum value of flow angle ratio product - >tp 
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25. Solve the equation for perforated plate (gauze) flow angle ratio - ‘'p 

2«p^ + - 1.21 Up ^ 2/3 
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26. Compute resistance of perforated plate (gauze) 

-2 


- K, 


Kp = 1.21 Hp 


- I 


27. Canpute flow angle ratio of hone.vcomb - '’n 


"p ‘"‘n 


‘n 


28. Compute resistance of honeycomb - K 

- 1 /4 

K|| = - ,201 Rt' Cn iXjj 
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steady Distortion Design 

29. Determine maximum azimuthal velocity in vortex - ^'2A 

30. Ccjupute maximum value of flow angle ratio product - ''‘n 

- 

"‘p "'ll 11 2 ;^ 

31. Determine maximum pre-contraction streamwise velocity deficit - 

32. Solve the equation for the perforated plate (gauze) flow angle ratio - 'p 


2a, 


+ \ ^ -1.21 V'^’aF 


+ 1.21 


IS 




33. 


<11 --n 'i,j 

1 I .201 j^ ;n a,, - ill '.pj 

Compute resistance of perforated plate (gauze) - K_, 




Kp - 1.21 .p 


_ 


34. Compute flow angle ratio of honeycomb - 

■p ‘11 

'M ' .7 

35. Compute resistance of honeycomb - 

K„ - .201 1 ;,.-' ... 

Conservative Design 

36. Canpare resistance of perforated plate (gauze) for atmospheric 
turbulence and steady distortion designs (Items ?6 and 33). Choose 
the larger - 

37. Canpare flow angle ratio of honeycoin:’ for atmospheric turbulence and 
steady distortion designs. (I tens 2 7 t.nd 34). Choose smaller - 'n 


38. Canputo honeycomb length to diameter ratio - ^ 


IbO 


39. Determine perforated plate open area ratio from Baines and Peterson 
pi ot F igure 114, - 

Perforated Plate Transmission Loss 

40. Determine the maximun mach number of flow incident on the ICS - 

41. Plot perforated plate thickness, I p, against hole diameter dp 
using the transmission loss criterion at the 24th one - third octave 
band. 


Ip ^ 5,39.10 
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5.04 


-. 1/2 
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dp meters 


42, Select a perforated plate thixkness and hole diameter consistent with 
the estimate of dj^s in I ten 10. 



Figure 114 Perforated Plate and Gauze Pressure Drop Coefficients 
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6.4.5 The Application 


of the ICS Design System to the JT9D 0T15D 


The step by step ICS design procedure is applied in Figures 115 and 116 for 
the JT9D and JT150 engines. A comparison between these designs and actual 


ICS's constructed for 

these engines is shown 

bel ow 


JT90 

Present 


Boeing 
3. 66m 

Nominal 
Radi us 

Desi gn 
3.47m 

P&WA 

3.66m 


Detail 

Dimension 

(Honeycomb) 

•0067m 

.0095m 

.0032m 

Detail 

Dimension 

,003m 

.0048m 

.001 57m 

(Perforated 
Plate or Gauze) 




Thi ckness 
(Honeycomb) 

.0117m 

,0762m 

.0381m 

Thickness 

.0009m 

.0031m 

.00079m 

(Perforated 
Plate or Gauze) 




Open Area 

54% 

51% 

46% 


Ratio (Perforated 
Plate or Gauze 


The most noticeable difference between the present design and the P&WA and 
Boeing ICS designs is the much larger honeycomb length to diameter ratios in 
these desigis - P&WA 8, Boeing 12, present 1.75. The P&WA and Boeing designs 
are very conservative in this respect. In other characteristics of the design 
the PSiWA and Boeing designs are generally conservative. 
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Figure 116 TCS Design Process for the JT15D 
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JT15D 


Nominal Radius 

Pres ent 
Design 

NASA 

Lewis (ICDl) 

Detail Dimension 
(Honeycomb) 

.00238m 

,0063m 

Detail Dimension 
(Perforated Plate 
or Gauze) 

.003m 

,0064m 

.0013m 

Thi ckness 
(Honeycomb) 

.0065m 

.05m 

Thickness 
(Perforated Plate 
or Gauze) 

,0009m 

.0006m 

Open Area Rati o 
(Perforated Plate 
or Gauze) 


m 

Again this comparison sliows a 
than would be required by the 

much greater honeycomb 
present design system - 

length to diameter ratio 
- NASA Lewis 8, present. 


2.31. The detailed dimensions of the screen material are also somewhat smaller 
for the present design than the NASA Lewis design. In addition, the open area 
ratio of the NASA Lewis gauze is more conservative than the present design. It 
should be noted that in both of these calculations estimates of the steady 
(quasi steady) distortiai field were made. Both vortical and axial velocity 
distortions were assimed to be of the order of the maximiri ambient wind speed. 

6.4.6 Summary of the Main ICS Design Characteristics 

0 The ICS should be located external to the engine in a region where, within 
the ambient wind test window, the engine induced flow is scmewhat greater 
than the maximua aiibient windspeed. This location is probably the best 
compromise between t'ne following requirements. 

a) ICS generated flow distortions to be small enough that no additional 
extraneous BPF tone energy is generated. 

b) Pressure drop across the ICS to be low enough so that the fan operat- 
ing point is unaffected. 
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c) Flow thi’ouqh the ICS to ho low enough so that acoustic transmission 
loss is negligible, 

cl) Steady or quasi-steady vortices to bo suppressed sufficiently far 

upstre<im of the inlet so that the residual axial velocity distortion 
can be suppressed by the flew contraction. 

The shape of the ICS is to seme degree flexible and while a basically 
hennsplieri cal desigi is suggested in this report, any shape that does not 
create inflow distortions is acceptable. 

Whether to mount the ICS on an impervious backplate or continue the screen 
to the engine nacelle has not been resolved. This question must be ad- 
dressed in conjunction with flight boundary layer simulation considera- 
tions . 

Discontinuities in the ICS should be minimized both for acoustic and flow 
distortion reasons. 

In view of the transfer of enerqy between transverse and streamwise velo- 
city conponents in a turbulent field subject to contraction, it is neces- 
sary to suppress all three velocity components. Consequently, the ICS 
should contain both a resistive elenent (perforated plate, gauze) and a 
torque withstanding elenent (honeycomb). If the inflow distortion field is 
dominated by steady or quasi-steady vortical elements and inflow turbu- 
lence is not a significant noise generator, the resistive elenent may be 
redundant, however in general this situation cannot be identified a priori. 
However, if this situation exists the use of a resistive element produces 
a conservative design. 

The relativi’ location of these two elanents is probably not important as 
regards distortion suppression. However, other considerations, such as 
pt’otecti on of the honeycomb, may suggest a preferred arrangement. 

The characteristi cs ot the ICS elanents are defined by the design system 
of Section 6.4.4. 


7.0 THE ASSESSMENT OF THE ICS DESIGN SYSTEM (TASK H) 


During the Boeing/Pratt & Whitney Aircraft Joint Noise Reduction Progran, data 
were gathered using pressure transducers mounted on the blades of a JT90 fan, 
with an ICS 1n place. These transducers were discussed in (29) and there it 
was indicated that it is not apparent, in general, which particular stimulus a 
BMT is responding to. This is especially true when an ICS is in place. In view 
of this situation the capacity for assessment of the ICS design systan using 
BMT data is reduced. It was however assuned that the BMT with an ICS in place 
was responding to the inflow velocity field and the assessment procedure was 
exercised. In addition, hot film data taken during the JNRP provided another 
and independent opportunity for assessment. 

7.1 The Assessment Procedure 

The assessment is made by cor paring a theoretical prediction, using elanents 
of the ICS design systen, with an estimate obtained from measured data. The 
quantity that was compared was the one dimensional PSD of the upwash velocity 
component in the azimuthal, k 2 * direction. A theoretical prediction of this 
PSD is directly obtainable from the ICS design systen. The estimate of this 
PSD, }, using BMT data requires sane discussion. 

The BMT samples the velocity and pressure fields along a spiral path, at an 
angle 0 to the fan face, in the mean flow fixed coordinate system. Figure 
117. It responds to duct fixed pressure and velocity fields and by averaging 
the signal on a once per revolution basis, this duct fixed or steady field may 
be determined. This and other operations were performed digitally. The steady 
contribution to the BMT signal is subtracted from the total signal to produce 
the duct unsteady field pressure signal. This duct unsteady signal is assimed 
to be due to perturbations in the incoming velocity field. The PSD of the un- 
steady pressure signal, P(f), may thus be determined. 
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iJNSTEAlJY 

PHissunr 


BMT SIGNAL SAMPUEO 
AT MOT FILM PROBC 
LOCATION (202” CCA) 



CmcUMFenENTIAL 
/ LOCATION OF 
BMT 


Figure 117 


Sampling of the BMT Pressure Signal 


No. as Shown in (28) the response function of the BMT has been estimated by 
using hot film and BW data gathered at the same circumferential location in 
the engine inlet. However, due to the sampling rate of the BMT, the upper fre- 
quency limit of the enpirical part of the response function was about 20 Hr 
and above this, a Sears function fonn was assmied. This is Justified, for in 
the frequency range of interest, the angle of Incidence of a given wave does 
not cause a significantly different response from the Sears function (see 
(45). It should be noted that while it is expected that the BW response func- 
tion Will have a similar form to the Sears function, the two functions are not 
directly comparable. The BW response function describes the surface pressure 
response while the Sears faction is a measure of the lift response of the 

olr,'. dPwn^h velocity 

PSD, i j^(f ) can be found. 
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The velocity of the blade at the BW radial location 
PSD,lj|(f), may be transformed into wavenimber space, 
lence assumption, so that 


is known and so the 
using the frozen turbu- 


1’ / k ^ ^ ~ 2j£ \ 

■’ \ \ ) 


M a 

l'-(f) 


n 


2n 
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Where k^ .ne wavenir,ber component along the spiral path a.d m is the 
BW flow relative Mach runber. Finally if the turbulence field is" cade up of 
mg eddies due to the contraction, the wavenunter directions are substan- 
tially aticuthal and consequently, ), can be determined 

' k 


O ( ''2 " coa,|. j " ’’n 


( ICj. ) cos<. 
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n this manner the required upwasb velocity PSO was obtained frm the measured 
m pressure si^al. The prediction of this PSD fr^ the ICS desi^ system 
enents was perfomed using equation 86. The input variables were those pre- 

at'sir ^ acquired. The computation was perfonned 

several values of the wave nunber component in the range of interest. 


The conparison between the theoretical prediction of the upwash velocity PSD 
F„ (kj) using Equation 105 and an estimate of the sane PSD usini BMT data 
as described above, is shewn in the tabulation below. 


•<2 



f^n (k?) 


m'l 

Es timate 

N/m 

Predi cted 


F rom 
BMT Data 


Using 

Equation 104 

2.31 
4.62 
6.93 
9. 24 
11.59 

1.56 
1. 71, 

3.45 
6.17 

8.45 

10“2 

10-2 

10-2 

10“2 

10-2 


1.6 10-^ 
5.43 10“5 
2.69 10-5 
1.52 10-5 
9. 79 1 0-5 
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The predicted values are very aiuch less than the values obtained free, the BMT 
data and the shapes are different. This is also true for estimates obtained 
from other BMTs. Thus in view of the reservations expressed previously about 
the BKT response this comparison is Inconclusive. 

In addition to blade pressure data, unsteady inlet velocity field infonoation 
was gathered from an array of split films mounted on a probe. This probe is 
described in (29). Both aaimuthal and streamwise velocity raeasuranents were 
made-inside and outside of the nacelle boundary layer with an ICS in place No 
infonoation was acquired spatially in the azimuthal direction and hence these 
data could not be used to estimate F„ (kj). However, the streamwise and 
azimuthal velocity component Intensities were found and are showp_in Figure 
118 (this data was previously published in (13)), from which an estimate of 
the upwash velocity cooponent intensity was obtained using Equation 23. Taking 
average values of the component intensities outside of the boundary layer the 
upwash velocity cooponent intensity was found to be approximately .5%. ’ 

A theoretical prediction of this quantity was made by using Equation 105 of 
the ICS design systao. A nonerical integration was performed over k, to give 
the theoretical value of the mean square value of the upwash velocity conpon- 
ent, froo which the intensity was computed and found to be approximately .25%. 

7.2 Conclusions 

0 Using BMT measuranents to determine the upwash velocity PSO. F (k ). 
there is a large discrepancy between this measured value and the vfilue 
predicted using ICS design systan elanents. This indicates errors in the 

design system and/or a source of excitation of the BMT in addition to the 
inflow velocity field. 

0 The agreement between predicted and measured values of the upwash turbu- 
lence intensity (obtained using hot films) is cmiparatively good. This 
provides support for the ICS design systan. 
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TURBULANCE INTENSITY (%> 



118 Component 
with ICS 


Turbulence Intenrsities Measured in 


the Inlet of a JT9D 


8.0 CONCLUSIONS 


A series of experiments have been conducted to evaluate the theoretical models 
to be used in the ICS design system. The major observations made during this 
test program were 

0 The flow contraction amplifies the transverse turbulent velocity com- 
ponent and at low contraction ratios, it suppresses the stream-wise 
velocity component. At high contraction ratios the streamwise turbu- 
lent velocity component increases, .alt hough it did not attain the 
starting value. 

0 Steady streamwi se -vel oci ty deficits are suppressed rapidly by the 
flow contraction. The azimuthal velocity component of a streamwise 
aligned vortex was substantially unchanged by the flow contraction. 
The streamwise velocity distortion associated with the vortex in- 
creased in amplitude and extent. 

0 The turbulence data gathered downstream of both honeycomb and perfor- 
ated plate was contaminated with self generated turbulence. The honey- 
comb screen suppressed the transverse turbulent velocity field. 

0 Steady streamwise velocity deficits are suppressed according to the 
resistance of the screen. Ln the resistance range studied here, the 
higher the resistance the greater the suppression. The azimuthal ve- 
locity field of the vortex was totally suppressed by the honeycombs 
and somewhat reduced and dispersed by the perforated plates. The 
streamwise velocity distortion associated with the vortex increased 
on convection through the honeycomb and the 66 percent open area per- 
forated plate. 
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0 The resistance of perforated plate Is adequately predicted by theory. 
The resistance of the honeyccmb is dependent on the entrance condi- 
tions and this precludes accurate estimates of the resistance for an 
arbitrary honeycaab. However an estimate may be made, by doubling 
Blasius' expression (Equation 6). 

0 The transmission loss of the thinnest honeycombs (.95 cm and 1.27 oa) 
is negligible and invariant with throughflcw speed. The thickest 
honeycomb (2.54 cm) has a transmission loss that increased with 
throughflokv speed to 2dB at 12.19 m/s. 

0 The transmission loss of the 66?S open area perforated plate is neg- 
ligible at the lower speeds (9.14 m/s and 6.1 m/s) but increases to 
2dB at 12.19 m/s. This attenuation is constant in 1/3 octave band. 

The 41% open area perforated plate transmission loss decreased with 
frequency and changes little with throughflow speed. The transmission 
loss is significant at the Icwer 1/3 octave bands. 

0 The presence of honeycomb panel corners and structure can disturb the 
radiation field of a simple discrete source by significant amounts. 
The higher the frequency the higher the disruption of the radiation 
f iel d. 

Subsequently the ICS design system models were modified to account for the 
results of the test program where necessary and these modified models were 
shown in Section 6. ^ transmission loss design criterion was established for 
the ICS and the effect of discontinuities in the structure on the radiation 
field was exanined. It was concluded that: 

0 Discontinuities in the ICS structure could affect the radiation field of 
discrete tones significantly but the broadband radiation field on a 1/3 
octave band basis would be substantially unnodified. 
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0 It was found that small source phase variations were sufficient to wash 
out the effect of the ICS discontinuities on the discrete tone radiation 
field, and it was suggested that controlled variations in engine speed of 
the order of 0.5 percent be made to accomplish this. 

The ICS design procedure is described on a step by step basis in Section 6.4 
and the procedure is exercised for the JT9D and JT150. The designs produced by 
this procedure are similar to those already constructed for those engines by 
Pratt & Whitney Aircraft, Boeing and NASA Lewis. The most significant differ- 
ence is in the honeycomb i/d ratio where the physical ICSs have a substantial 
overkill. The assessment of the procedure was made in Section 7 using BMT and 
hot film data. It was concluded that: 

0 The hot film data provided some support for the validity of the ICS design 
systen. 

0 The BMT data provided little support for the validity of the ICS design 
system. It was suspected however that the BMTs were responding to excita- 
tions in addition to the inflow velocity field which for the purpose of 
ICS design assessment were extraneous. 

The similarity between the ICS design produced by this procedure and actual 
ICSs is also encouraging in view of the proven effectiveness of those designs. 
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appendix a 


FLOW MEASUREMENT INSTRUMENTATION AND ACCURACY 


Measurement Type 1 . (Mean Velocity In the Duct) 

. In the duct (30.49 to 121.95 m/s), a wedge type 

For mean velocity measurements in the duct 

5 hole oressure probe was used. The probe has a 15 included angle. .-9 
width rt the trailing edge and .051 O" diameter pressure tap. The pro e as 
one total tap (P^) and two static taps (P, » P^ located one on each 
side). Two certified pressure transducers, 2.54 and 101.6 cm H 2 O F. ., 
used to cover the full range of pressures. 


A cmnplete calibration included velocities of 34.45. 78.96 « 

each velocity, yaw angles of t.5* and .!». This 

whole systan which includes the probe, transducer and voltmeter ‘Oe sam 

systan components were always used together throughout the tests). 

is shown in Figure 119. 



Fiqare U9 Wedge Probe Sg.rtem 
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The usual mode of operation of this type of probe is to rotate the probe until 
it is aligned with the flow stream (as determined by the two static tap 
readings being the same). Because of the small flow angle range expected in 
this test, the probe was operated in the fixed angle mode. The velocity 

(magnitude and angle) was determined from pressure readings (pTn-P« and 
P -P 1 ^ ^ 

Measurenent Type 2, (Turbulence in the Duct) 

For turbulence measuranents inside of the duct (24.39 to 121.95 m/s), X-array 
probes with hot-wire sensors were used. The probe tested was TSI model 1227 
with a 9 urn diameter, .279 cm long W wire sensor_The sensor was operated at 
300°F overheat using a TSI model 1053A anemometer. 

The systOT accuracy for X-array probes when used with linearizers and 
sum/difference amplifiers (to obtain turbulence components) is estimated to be 
+6.3%. The system is shown in Figure 120. 




X ARAAY HOT WIRE PROBE 
BpM TUNGSTEN WIRE 


ANEMOMETER 


L1NEARI2ER 



ANEMOMETER 



magnetic TAPE RfCORDtR 


Figure 120 Hot Wire Si)stem 


175 










Measuranent Typo 3. (Moan Velocity Ahead of the Duct) 

For moan velocUy measuranents ahead of the duct (2.44 to "^*J’,Firw1th 
elaoent hot-wire sensor was used. The probe tested was a TSl ^ 

, , « dlaneter. .569 ca, ion, Vi wire sensor. The sensor was operated 

300''f overheat using a TSi »del 1063A anemoneter. 

TO detemiine the accuracy, three pre-test and one post test 
erfon^ed. The caiibration velocity range was 1.83 tn 20.12 m/s and t r 
probes were involved. Based on 12 sanple cmpansons. the accuracy wa 
to be V2 8*. The post test and pre-test calibrations agreed wUh,n 2.n. 

applies to the probe. an..,eter and The syst. 

accuracy for X-array probes when used witt. linearizers and sun/differe^e 
^plifiers for mean velocity (magnitude a,d imgle) measurohents is tS.O* and 

+ 3 °. The systan is shown in Figure 121. 


SINGLE ELEMENT HOT WIRE PROBE 
DIA. WIRE 


o 


anemometer 


DC 

VOLTMETER 


F 1 


|J,> Hot M-iro Ftic-trw 


Measurement Type 4. (Turbulence Ahead of the Duct) 

For turbulence measuranents ahead of the duct (2.44 to 12.2 
probes with hot-wire sensors -ere used. The probe tested was a TSl «di 1-40 
with a 9 un. dUneter. ,110” long W wire sensor. The sensor was oper ... 
jCO-’f overheat using a TSl model 1053A ananoneter. 






The system accuracy for X> array probes when used with linearizers and 

sun/ difference amplifiers (to obtain turbulence components) is estimated to be 

+6.3%, The system is shown in Figure 120. 

A summary of measurement accuracies is given in Table AI. 


TABLE A-I 


Measurement Type Instruments Accuracy (2 s.d.) 


1. 

Mean Velocity in the Duct 

Wedge pressure probe, transducer. 

_5X, 

10 

2. 

. Tur bul ence i n. the Duct - 

- Same as type 4. Correlation factor 
of 1.36 must-be applied for scale 
effect. 

Same as 

type 4 

3. 

Mean Vjelocity Ahead of 
the Duct 

9pm X .220" Hot-wire probe, 
anemometer, voltmeter 

2.8* 




L inear i zer 

2* 




Sum/Difference Amplifier 

1* 




System; 

5.0*, 

30 

4. 

Turbulence Ahead of the 
Duct 

9 m X .110" Hot-wire probe, _ 
anemometer, voltmeter 

3.9** 




Linear! zer 

2*- 




Sum/0 if fer ence Amplifier 

1* 




System: 

6.3* 



*Estimate based on steady state calibration data and user experience 

Post test calibration agre> with pre test within l.TIt 
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APPENDIX B 


THE INVARIANCE OF THE POTENTIAL VORTEX 
VELOCITY FIELD UNDER AXISYMMETRIC CONTRACTION 

Applying the conservation of circulation to particles in a vortex subject to a 
circular axi synunetri c contraction (Figure 122) produces 



Figure 132 Thf Axisuimetric Contraction of a Vortex 




as 


A 





as 


B 


( 'a ) 


A 


a 




(’u ) 


B 


119 


120 


but for a potential vortex circulation is constant at any radius 


^'.>u ( ' lO ' ’’.’a ( ‘ a ) ''a 
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Thus substituting for | ^ in 120 


l 'a I ( 'a ) 


12 ? 


i.e., the azimuthal velocity field is invariant. 
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APPENDIX C 


THE MODIFICATION OF THE VORTEX STREAMWISE VELOCITY 
ON CONVECTION THROUGH HONEYCOMB AND PERFORATED PLATE 


a) Honeycomb 


Applying Bernoulli's equation across the airfoil (following Batchelor) 
yields (Figure 123): 



AIRFOIL MONIYCOMB 

OR 

PERFORATrO 

KAU 


Figure 123 The Convection of a Vortex Through Honeycomb 

If the radial velocity component is neglected then 


u 






124 


throughout the vortex, ah is the total pressure loss due to viscous 
effects at the airfoil surface. The static pressure distribution in the 
vortex, , is a function of the azimuthal velocity field. Equation C2 
illustrates the point made in the text that the axial velocity field 
depends on the relative magnitudes of the pressure loss across the 
airfoil and the azimuthal velocity field. 
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2 2 
IJ, . - II 
1 ‘A 


Then If = u - » Equation 124 becomes 


r ( '’■■■ ■ i’a ) " '^' 2 ^ * ;t -'h 


and if ^ == ’J.. 


\u - ^ I 2 , . 

^"lA - 2iT I - r ( - !=a ) + 


i 

‘ '■ I 


Considering now the flow across both airfoil and honeyccmb (see Figure 122) 
and again applying Bernoulli's equation gives 


P„ ♦ 5 •• U/ - AH - All^ - l>B + ^ OjB* 

where Is the pressure loss through the honeycomb. 

Solving for , 


!■» I • - r ( ‘ 'Ii.s ) 


If again the deficit Is assumed small 


I I 2 


( " f’u) ' 7 ( ^ ) j 
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Now If Equation 127 is subtracted fron 130, 


AUj„ - ,UJ, 


A ' ri“ I ;■ ( i'a - "I, ) - 


2 2 
. ~ AIlj, 


131 


or 


AU 


in 


AU 


1 + - 


u. 


lA 


( u,, AU*“ j t’n " f’A + 


} "'2A 

S j ■ 2ir~MJ 


lA 
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Assuming that the total pressure loss across the honeycomb is constant and 
equal to the loss at . 


Then ) 


133 


Downstream of the honeycomb, the azimuthal velocity component is zero and 
consequently the static pressure gradient cannot be maintained and the static 
pressure will become uniform so that 

^ 134 

Substituting Equations 133 and 134 into Equation 132 gives 


AU 


IB 


AU, 


= 1 + 


U, 


'lA 

Now in a vortex 


Au^^ I Pa- ■ Pa 1 


P; 


'2A 
7U AU 


lA 


dp ^ ''^2A 
dr 


so that 


t3 


A' 


->a - - ... 
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136 


137 
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Substituting this expression in equation 13b 


AU 

Afi; 


IM 

A 


n AIL 


A 


^ f - "if' I 
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( r • ^ ( 

The integral is always positive, the initii.l deficit u, s 

aencit AUj^ may be positive 

(wake) or negative (jet). The honeycomb thus amolifip. . . 

and suppresses or Inverts a vortex Jet In hnth deficit 

changed. distortion shape is 


For the maximim axial deficit (excess) rati 

(the .easored aai^tha! veiocuT ^“idTcaT:: 
a Kankine Vortex) so that 


0 across the honeycomb, assume a 

approximated by 


/ 


u 


2a 


cir u- 


o 2 a in.ix- 

and the ratio becomes 

AU 
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AU 


IB ma x 
lA max 


tL, 


2 A max 


A max 
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d) Perforated Plate 


If the azimuthal velocity field i 


perforated plate Equation 120 is replaced by 


s not destroyed as in the case of the 


K. > 7 i ti 




and so for a small deficit, 


3 ■ (“ll/ " "2,,' ) 


AU, _1 ( 

■ ir / '» P.. • nil , I 


U,., ' 
2U 
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is: 


14? 


Recalling Equation \?.l and subtracting it from Equation 142* 


now 


''“in • ''“lA ’■ rV ! I’n - Pa ■* ''"k I + " 


^2B ’ 


(”2B^ - “2A^) 


Thus 


4U 


IB 


au 


= 1 + 
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2 / 2 


lA 




pp “ p. + AH 


2A 


(„^-i) 


s 1 -TD mT 


lA 


Again assuming AHg == ( ~ ) 

AUi 


'IB 


AU 


= 1 + 


lA 


U 


2A 


P U AU,- 
lA 


(n^-l) 


I ( ^Ao' " Pa ) ‘ ( “ Pb ) 


2U AU 


2A 


If further a Rankine vortex is assumed on both sides of the perforated 
then 


^^^IB max 

, . “2A^ma. (l ' 

^^lA max 

''“lA „,ax 


143 

144 

145 

I 

1 

146 

plate, 

147 
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APPENDIX D 


THE UPWASH VELOCITY ELEMENT OF THE SPECTRUM TENSOR AND THE TARGET (FLIGHT) ONE 
DIMENSIONAL PSD AT THE FAN FACE. 

The upwash velocity, assuming a frozen turbulence field, may be determined 
from a simple rotation of the Cartesian coordinates. Figure 124, and written as 

u (x) = - u, (x) sin (x) cos 


Figure 124 The Turbulence Upwash Velocity Component 


If the three dimensional correlation is computed and Fourier transformed, the 
relationship between the spectrum tensors is found to be 


(k) “ (k) sin' 




2 14 

(k) cos . ^ + cross terms 






Now the one sided one dimensional PSD over the ranges 

' ’^ 21 , " *^2 ' '^211 * " ^311 



Is 
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The integration is first performed by approximating the integrand with 

its value at the mid point in the wave nunber range, so that 
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substituting this in Equation 152 yields 

C'„(k2)^32N 
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/ 


-k, 


4k2^sin^'' .1- 


^ ( ^3u ■'■ cos^ Yq 


But over the integration range 


4 



+ 



■> s 
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and 


4k, 


sin Y 



N N 



COS 
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The integral in Equation I 54 


can then be approximated so that 



in the ranges 

is then the target 
with an ICS. 


* 

^ ■ '‘21, - kj ■: kjy , kj . 

(flight) PSD that fs to be aimed for 


*^3 ' ''3U Eguation 158 
on the static test stand 
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